Two-mode back-action-evading measurements in cavity optomechanics 
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We study theoretically a three-mode optomechanical system where two mechanical oscillators are 
coupled to a single cavity mode. By using two-tone (i.e. amplitude-modulated) driving of the cavity, 
it is possible to couple the cavity to a single collective quadrature of the mechanical oscillators. In 
such a way, a back-action-evading measurement of the collective mechanical quadrature is possible. 
We discuss how this can allow one to measure both quadratures of a mechanical force beyond the full 
quantum limit, paying close attention to the role of dissipation and experimental imperfections. We 
also describe how this scheme allows one to generate steady-state mechanical entanglement; namely, 
one can conditionally prepare an entangled, two-mode squeezed mechanical state. This entanglement 
can be verified directly from the measurement record by applying a generalized version of Duan's 
inequality; we also discuss how feedback can be used to produce unconditional entanglement. 

PACS numbers: 03.67.Bg,42.50.Lc,85.85.-|-j 



I. INTRODUCTION 

The idea of a back-action-evading quantum measure- 
ment is by now a well-known concept [iHl] . The simplest 
and most studied realization is the continuous position 
measurement of a simple harmonic oscillator. If one tries 
to simultaneously follow both quadratures of the oscil- 
lator's motion (i.e. both the amplitudes of the sine and 
cosine components of the motion) , then the unavoidable 
effects of quantum back-action imply that one cannot 
improve the measurement indefinitely by increasing the 
measurement strength. This leads to the quantum limit 
on continuous position detection [31 13]. If instead one 
measures just a single quadrature of the motion (say the 
X quadrature), there is no such limit as one has "evaded" 
the back-action: the measurement back-action only heats 
the conjugate P quadrature, which has no effect on the 
measurement as the dynamics does not couple it to A". 

Such back-action-evading (BAE) measurement 
schemes allow one to measure one quadrature of a 
narrow-bandwidth force acting on a mechanical os- 
cillator with arbitrary precision. They also naturally 
lead to squeezing of the measured quadrature, albeit 
a conditional squeezing, where one must have access 
to the full measurement record to produce uncondi- 
tional squeezing. Back-action-evading schemes are also 
known as "quantum-non-demolition (QND) in time" 
measurements, as one is making a QND measurement 
of an observable which is explicitly time-dependent. 
These ideas have recently been discussed [5] within the 
specific context of quantum optomechanics [7], where 
a mechanical oscillator is dispersively coupled to a 
driven cavity. A back-action-evading single-quadrature 
measurement was even implemented experimentally 
by Herzberg et al., using a nanomechanical oscillator 
coupled to a driven microwave- frequency cavity fS" . 

One could now ask about exploiting the idea of back- 
action-evasion in more general settings. The simplest 
generalization would involve two harmonic oscillators. 



implying four quadratures of motion. Such a two-mode 
BAE scheme would involve measuring two quadratures 
of motion, with the back-action only driving the unmea- 
sured conjugate quadratures. If one could suitably couple 
the system, one could imagine measuring both quadra- 
tures of an applied force without any quantum limit. 
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FIG. 1: (Color online), (a) The system studied consists of 
two mechanical oscillators, each coupled to a cavity (or cir- 
cuit) mode. Measurement is performed via phase-sensitive 
detection near the cavity resonance of the cavity output field. 
(b) Frequencies in this system, defined with respect to the 
cavity resonance frequency uic- The blue lines indicate the 
standard mechanical sidebands, at ±uia and ±0^5. The cavity 
is driven at ±u}m ~ ±(a;a -|-aji,)/2, as indicated by vertical red 
lines on either side of the cavity resonance, (c) Representa- 
tion of the two-mode BAE scheme in the ideal, fully symmet- 
ric case. The observable X+ is directly measured, and the 
observable P+ is heated by the corresponding back-action. 
P- is dynamically coupled to X+, and therefore effectively 
measured provided that f2 2> 7. In this limit, X- (which is 
dynamically coupled to P_ ) is also heated by the back-action 
of the measurement. 



That such a scheme is possible was recently discussed 
in detail by Tsang and Caves O [10] , though similar ideas 
have been discussed earlier (see, for example, Ref . [TT] or 
App. D of Ref. T^. From a practical perspective, it ap- 
pears challenging, as it requires one of the two oscillators 
involved to have effectively a negative mass (and hence a 
negative frequency). Similar to standard single-oscillator 
back-action-evasion, one might expect that there is also 
a squeezing aspect to such two-mode BAE. This was dis- 
cussed by Hammerer et al. [13] : the squeezing here of the 
delocalized measured quadratures naturally gives rise to 
an EPR-style entanglement between the two modes. 

In this paper, we now extend the above idea of two- 
mode back-action-evasion to quantum optomechanics, 
specifically considering a system where two mechanical 
oscillators are coupled to a single electromagnetic mode, 
as depicted in Fig.[T] Such two- mode optomechanical sys- 
tems have recently been realized experimentally [14 . We 
show that the two-mode BAE scheme can be realized in a 
relatively simple fashion in such a system, using a gener- 
alization of the two-tone driving scheme used to achieve 
single- mode back-action-evasion [IjiGj. This is in contrast 
to previous proposals, where to achieve the requisite neg- 
ative mass, one of the two modes was not mechanical, but 
rather an atomic ensemble |13j or another cavity mode 
[5]. Unlike the single- mode back-action-evasion case, we 
expect that this scheme will not be susceptible to para- 
metric instabilities [TS] , since here the number of photons 
in the cavity will not oscillate at double the mechanical 
resonance frequencies. 

Note that applications of BAE measurement ideas to 
a system of two mechanical oscillators have previously 
been studied [TBH^ . motivated by attempts at gravita- 
tional wave detection. The approaches presented in those 
works differ significantly from our scheme; in particular, 
they involve an explicit coupling between the mechani- 
cal oscillators, something that is not required with our 
scheme. 

We provide a thorough discussion of the optomechan- 
ical two-mode BAE scheme, including a discussion of 
expected experimental imperfections, and methods for 
countering these. In addition, we provide a full analysis 
of the force sensitivity of this scheme and the possibility 
of beating conventional quantum limits. Unlike previous 
work [9l[T0], we include dissipation in this analysis. 

Finally, we also consider the generation of station- 
ary, conditional entanglement using this scheme, and the 
possibility of turning this into unconditional entangle- 
ment using feedback. This is in contrast to more con- 
ventional measurement-based entanglement generation 
schemes, which rely on strong measurements and post- 
selection; such schemes are well-established in optics [23|- 
[25] , and have recently been discussed in the context of 
superconducting circuits [261 El] ■ Our scheme also differs 
from that in Ref.[T3] as we generate stationary entangle- 
ment, and do not require a single-shot strong feedback 
operation. 

The remainder of this paper is structured as follows. 



In Sec. [TT] we derive the Hamiltonian and Heisenberg- 
Langevin equations describing this system and introduce 
a method for compensating for system asymmetries. In 
Sec. |III| we calculate the back-action heating effect and 
the measured noise spectrum, with a view to assessing 
the force sensing limits of our system. 

In Sec. |IV| we rigorously compare the force sensitiv- 
ity of our scheme, including imperfections, against con- 
ventional quantum limits. The situation is significantly 
more complicated when one includes the presence of dis- 
sipation (unlike Refs.lUland fTU)) . We find that for a signal 
force applied at the mechanical resonance frequency, the 
number of added noise quanta can be made to go to zero 
(below the quantum limit of 1/2), 
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even in the presence of dissipation. However, this is only 
possible in a narrow bandwidth about the mechanical 
resonance and requires a careful tuning of optomechan- 
ical couplings. For a signal force off-resonant from the 
mechanical resonance frequency (detuned by an amount 
A), we find that the added noise can be made to go to 
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well below the standard quantum limit in this case 
(^add[wQ -I- A] = A/7 where 7 is the average mechani- 
cal damping rate) and even below the full quantum limit 
(i.e. allowing for optimal detector noise correlations). 

The conditional dynamics of the mechanical oscilla- 
tors, under a continuous measurement via the coupled 
cavity, are described in Sec. [V] The steady-state con- 
ditional variances are determined, and these are used to 
determine in which regimes the mechanical oscillators are 
conditionally entangled. The presence of all-mechanical 
entanglement is determined from Duan's inequality [28j . 
which here (in the case of symmetric optomechanical cou- 
plings and mechanical damping rates) becomes the re- 
quirement 



Vx 



Vp_ = 
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where Vx^ {Vp_ ) is the variance in the sum (difference) 
of the mechanical X (P) quadratures, nth is the effec- 
tive thermal occupation of the mechanical environment 
(see Eq. (83al), rj is the quantum efficiency of the mea- 
surement (see Eq. ([73])), and C is the so-called optome- 
chanical cooperativity, essentially quantifying the mea- 
surement strength (see Eq. (29l). Both quantities on the 
left-hand-side of Eq. (pj) may be obtained in a straightfor- 
ward manner from the measurement record. Therefore, 
the mechanical oscillators should be entangled provided 
that the cooperativity obeys 



C> 



2(nth + 1/2) 
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From an experimental perspective, this would appear to 
be a modest requirement, well within the reach of cur- 
rent optomechanics experiments based on microwave cir- 
cuits [25^. In Sec. VI the total unconditional variances of 
the collective quadratures are calculated, and a feedback 
scheme is introduced that allows the conditional entan- 
glement to be converted to unconditional entanglement. 



II. SYSTEM 
A. Hamiltonian 

The system (see Fig. IT]) is composed of two mechanical 
oscillators, with resonance frequencies uja and Uh, each 
independently and dispersively coupled (with strengths 
Qa and (7b, respectively) to a common cavity mode having 
resonance frequency ujc- The Hamiltonian is 

T-l = Wafl'a + ajb& & + '^cC c + g^j (a + a* ) c' c 

-h56(& + 6^)c^C-Hi?drivc+i^diss, (5) 

where a and 6 denote mechanical mode lowering oper- 
ators, c denotes the electromagnetic mode lowering op- 
erator, and Hdrive accounts for driving of the electro- 
magnetic mode. The term iJdiss accounts for dissipation, 
with the modes subject to damping at rates 7^, 7^ and 
K, respectively. 

To realize the two-mode back-action-evading dynam- 
ics suggested by Tsang and Caves |9j, we need the two 
mechanical oscillators to have equal and opposite fre- 
quencies, and moreover, have the cavity (the detector 
in this scheme) only couple to a single collective quadra- 
ture of the two mechanical oscillators (an operator which 
is explicitly time-dependent). Both these requirements 
can be realized by simply adapting the two-tone (i.e. 
ampltiude-modulated) driving scheme used for ordinary 
single-quadrature BAE in Refs. IHI and |H1 One applies a 
drive to the cavity at both Wc + w™ and Wc — tOm, where 
Wm = {ua + i^b)/'^ is the average of the two mechanical 
frequencies, i.e.: 



[a = a, b,c): 



H, 



drive 



= (^;^ 



-K£:*e*""*)e*"'^*c-Kh.c. (6) 



Working in an interaction picture with respect to Hq = 

uj,yi ( a^a + b^bj +lUcC^c, one finds that the frequencies of 

the mechanical oscillators are now ±D, as desired, where 
2fl — uja — ^b- We further specialize to an optomechani- 
cal system in the good cavity limit (a;„j 3> k) and assume 
large driving amplitudes \£+\, \£-\, such that one can lin- 
earize the optomechanical interaction. In this standard 
regime, one finds (as desired) that in the interaction pic- 
ture, the cavity only couples to a collective mechanical 
quadrature with a time-independent coupling. Specif- 
ically, defining the the quadrature operators for each 
mode in the interaction picture in the standard manner 



Pa = -i(a-a^)/V2, 



(7a) 
(7b) 



X± — (Aa±A(,)/v2, 

P± = {Pa±Pb)/V2, 



and collective (canonically conjugate) mechanical 
quadrature observables via 



(8a) 
(8b) 



n = n IX+X^ + P+P-] + GX+Xc - GaX^X^ 

+Hdiss- (9) 

Here, the many-photon optomechanical couplings are 
given by 



the optomechanical Hamiltonian takes the form 



G = V2{ga+gb)c, 
Gd = \/2(gb - ga)c, 



(10a) 
(10b) 



where c — \E±\/ ^lu"^^ + k^ /4, is the steady-state ampli- 
tude of the cavity field at the driven sidebands. Note 
that we want to be in a regime where k > fi, such that 
the cavity responds fast enough to monitor the collective 
mechanical oscillation. The details of this derivation are 
given in App.|X] We stress here that the interaction terms 
in Eq. ^ , even accounting for contributions from the off- 
resonant sidebands of the driving fields, will not include 
terms oscillating at double the mechanical resonance fre- 
quencies. Accordingly, the cavity photon number will not 
oscillate at this frequency and the parametric instabili- 
ties associated with single-mode back-action-evasion [T3] 
should be avoided in the two-mode case. 

For equal optomechanical couplings (i.e. Gd — 0) and 
equal mechanical damping rates (7a = 7b = 7) this sys- 
tem perfectly realizes the BAE scheme of Tsang and 
Caves. To see this, consider the Heisenberg-Langevin 
equations for our system in this ideal symmetric limit. 
Defining the vector of collective mechanical quadrature 
operators as 



F=(X+,P_,A_,P+ 



a standard calculation yields 
d 



dt 



V ^M-V + Fi 



BA 



N-e 



(11) 



(12) 



Here, M describes the oscillator dynamics in the rotating 
frame, 



M 



-7/2 r2 
-n -7/2 
-7/2 n 
-r^ -7/2 



(13) 



while Fba is the back-action force from the cavity: 

FBA = (0,0,0,-GX,f. (14) 

Further, 

^= (X+an, P-.in, ^-an, P+Anf, N = ^U, (15) 

describe the Langevin noise arising from the dissipative 
baths of each mechanical resonator, with I„ denoting the 
n X n identity matrix. Finally, the equations of motion 
for the cavity quadratures are: 



2 



^/kXc 



Pr 



-GX+ - -Pc + V^PcA 



(16a) 
(16b) 



We thus have the desired behaviour, namely: 



• The collective X^ and P_ quadratures are linked 
dynamically the same way as X and P for a har- 
monic oscillator, but are dynamically independent 
from their conjugate pair (P+ and X_). In the 
language of Tsang and Caves, the commuting ob- 
servables X^ and P_ form a quantum-mechanics- 
free-subspace. 



• From Eq. (16b), we see that the cavity mea- 
sures X+ (and hence the dynamically linked P- 
quadrature), and its back-action only drives the 
unmeasured mechanical quadratures P+ and X_ 



(C.f. Eq. (14)) 



One thus has a setup which allows the back-action- 
evading measurement of both the collective quadratures 
X+ and P_ . Information about the motion of these col- 
lective quadratures appears at the cavity output at side- 
bands ±J7 from the cavity resonance. 



B. Compensating for unequal optomechanical 
couplings 

While the scheme is easy to understand in the perfectly 
symmetric case, any real expeirment will have to contend 
with both an asymmetry in the optomechanical couplings 
(ga 7^ gb) and an asymmetry in the mechanical damping 
(7a 7^ 7b)- Asymmetric couplings cause Gd in Eq. ([9]) 
to be non-zero, implying that the cavity measures both 
the X+ and X^ collective quadratures. As a result, the 
measurement is no longer back-action-evadig (i.e. the un- 
wanted measurement of X- causes back-action to drive 
P_, which then corrupts the measurement of X^). Even 
if there is no coupling asymmetry, if 7^ ^ 7f,, the damping 
terms in the Heisenberg-Langevin equations will link X^ 
and X- , also ruining the back-action-evasion. Both these 
effects can easily be seen by writing out the Heisenberg- 
Langevin equations in the general asymmetric case; see 
App.lBl especially Eqs. ^B2\ and ([B5|. 



We will show in later sections that while small devi- 
ations from the ideal symmetric system cause a depar- 
ture from perfect back-action-evasion, one can still beat 
quantum limits on detection as well as generate squeezing 
and entanglement. Here, we point out that it is possible 
to exactly null the deleterious effects of asymmetric op- 
tomechanical couplings by introducing additional para- 
metric drives on each mechanical resonator. Note that 
attempting to compensate for the asymmetric optome- 
chanical couplings by modifying the weights of the two 
cavity drives only leads to more unwanted interaction 
terms. 

To begin our analysis of the case ga ^ 5fc, it is useful 
to introduce a new set of canonically-conjugate collec- 
tive mechanical quadratures, such that in this basis, the 
cavity only couples to a single (collective) quadrature. 
Defining the rotation angle 9 as: 



= tan-i(G'd/G), 



(17) 



we introduce: 






x± 


= cos 6 X± =F sin 9 X^ , 


(18a) 


P± 


= cos6l P± =Fsin6' P^,. 


(18b) 



The Hamiltonian of Eq. ^ , in terms of the observables 
introduced in Eqs. (18a[) and (|18b|, takes the form 



n = n[^x+x^ 

where we have defined: 



Xi 



P+P-]+ GX+Xe 



Pl-Xl-Pll, 



H. 



diss 



p 



tan 26* 



2Gd/G 
1 - Gl/G^ ^ 2gagb 



9 9 

9a - at 



n = r)cos20 = r!^-9/£ = 2 ^"^^^ 



l + Gl/G^ 



.9a + d ' 



G 



Gsec6i = GJl + G2/G2. 



(19) 

(20a) 
(20b) 
(20c) 



One can recast the Heisenberg-Langevin equations de- 
scribing the system in terms of these rotated observables. 
The system retains the form of Eq. ( 12 1, now in terms of 
the vector of rotated observables. 



F = (1+,P_,X_,P+)^, 
with the input noise vector given by 

C — (^+,m5 -P-,irn A"- jn, P+Jn) ■ 



(21) 



(22) 



The appropriate matrices are specified in Eqs. (B9) and 

dBTol). 



The cavity is now only coupled to a single collective 
quadrature A+, with a modified coupling constant G. 
The free evolution of the oscillators in this quadrature 
basis is, however, not of the desired form; the back- 
action-evasion is ruined by the terms in the last line 



of Eq. (19), which dynamicahy Hnk the measured ob- 
servables (X+,P_) to their (perturbed) conjugate pair 

It is now natural to ask whether it is possible to mod- 
ify the mechanical system in some simple way so as to 
eliminate the last line in Eq. (19). Such a modification 



would require introducing both position and momentum 
couplings between the two mechanical oscillators, some- 
thing that is not experimentally feasible. A more mod- 
est approach would involve simply adding terms to the 
Hamiltonian of each mechanical oscillator (but not intro- 
duce any coupling). We consider parametrically driving 
each mechanical oscillator, and thus adding a term /fpa 
to the system Hamiltonian of the form: 



i?pa - -A ( a2 + at2 + 52 ^ ^t2 



(23) 



This is written in the interaction picture; in the 
non-rotating laboratory frame, this corresponds to an 
equal-strength parametric modulation of each oscillator's 
spring constant at the frequency 2a;,„. Writing this in the 
tilde quadrature basis of interest, we have: 



H„ 



A P' 



PI 



xi-X' 



(24) 



We thus see that by taking A = prj/2, these terms com- 



bine with the last term in Eq. ( 19 ) to give a Hamiltonian 



which has the general back-action-evading form: 



nc = n + H^^ = n{ x+x_ + p+p_ + gx+x, 



^ + ^^c 



pn(pl- X] 



H, 



diss- 



(25) 



Note that the required value of the parametric mod- 
ulation corresponds to requiring the amplitude of the 
spring constant modulation of each oscillator to satisfy 
Aki/ki = Apfl/tJi (i,= a,b). In this Hamiltonian, the 
measured observables X^ and P_ are dynamically inde- 
pendent of their conjugate quadratures as required. The 
extra terms in the last line of Eq. ( 25 ) only modify the dy- 
namics of the unmeasured subsystem. We thus see that 
two-mode BAE is possible even without having perfectly 
symmetric optomechanical couplings. 



III. NOISE SPECTRA 

To gain a better appreciation of the two-mode BAE 
scheme introduced above, we calculate the noise spec- 
tra of the measured and perturbed collective mechanical 
quadratures, and of the homodyne current from measure- 
ment of the cavity output. These may be calculated from 
the Heisenberg-Langevin equations in the usual man- 
ner ^; the details are given in App. [C] Note that all 
quantum noise spectra quoted throughout this paper are 
symmetrized quantum noise spectra. Of particular in- 
terest here will be the effects of asymmetries (either in 



the optomechanical couplings or in the mechanical damp- 
ing rates), which will cause deviations from the perfectly 
back-action-evading result. Readers not interested in the 
effects of imperfections (i.e. asymmetries) may wish to 
read Sec. Ill A and then jump to Sec. IV where we dis- 
cuss beating the quantum limit on force sensing. 



A. Ideal symmetric case 

We start with the ideal case of symmetric optome- 
chanical couplings {ga = gb) and symmetric mechanical 
damping (7^ = 7b = 7). The symmetrized noise spectral 
density of an observable Z = [V]i is defined by 



•SzM 



+ 00 



dte' 



^* {{Z{t) JiO)}). (26) 



For the measured observable X+ it is: 



5'x+H = ^{na + nb + 1) 



7/2 
(7/2)2 + (^ _ ^)2 



7/2 
(7/2)2 + (c^ + r!)2 



(27) 



where n^ and nj, are the temperatures of the two mechan- 
ical baths (expressed as a number of thermal quanta). 
We are working here in the interaction picture as always 
(i.e. uj is measured with respect to the average mechan- 
ical frequency Wm), and have assumed the good cavity 
(or "resolved-sideband") limit k <^ Wm. The spectrum 
of Eq. (27) consists of two Lorentzians, centred at ±fl, 



and weighted by the thermal and quantum fluctuations 
of the mechanical oscillators. Since we are in the perfect 
back-action-evasion limit, there is no dependence on the 
optomechanical coupling G. Turning to the perturbed 
observable P+, the spectrum is given, in the limit k ^ Q, 
by: 



Sp+ [w] 



1 



la + nb + l + C{2nc + 1) 



7/2 
(7/2)2 + (^ _ ^)2 

, 7/2 

' (7/2)2 + (c. + 0)2 



(28) 



where fie is the number of thermal quanta associated with 
cavity dissipation, and we have quantified the optome- 
chanical coupling G in terms of the cooperativity param- 
eter [30], defined as 



C^ 



2G2 

7K 



(29) 



The cooperativity here is essentially the measurement 
strength, and thus describes the back-action heating rate 
of the perturbed collective quadrature. It is equiva- 
lent to the quantity tt-ba used in Ref. [B] to denote back- 
action heating of the perturbed quadrature (as a number 



of quanta) in optomechanical single-mode back-action- 
evasion. 

We now turn to the spectrum of the cavity output. 
We consider the case of a single-sided cavity, with homo- 
dyne detection of the cavity output field. The measured 
spectrum of X^ is then defined in terms of the total ho- 
modyne current noise spectrum Sj[uj] as: 



cnicas 



[lj]^Si[u]/IC^ 



(30) 



where JC = 2_BG/\/k/2 is a measurement gain with B 
denoting the amplitude of the local oscillator in the ho- 
modyne detection scheme. Further details of the calcula- 
tion of this quantity are given in App.O Taking the most 
interesting limit k > 17,7 (i-^- the cavity is fast enough 
to measure the collective quadrature), we find 



S 



x^ 



M 



Sx+ [w] 



+ ^(2ne + l). 



(31) 



The second term here is the imprecision noise of the mea- 
surement (i.e. output shot noise). In the perfect back- 
action-evading limit considered here, it can be made ar- 
bitrarily small without any resulting back-action penalty 
by increasing the cooperativity via the driving strength 
£±. Note that the standard quantum limit result at res- 
onance (applicable when one is limited by back-action) 
is that the added noise must be at least as large as the 
zero-point motion of the oscillator. We thus see that one 
beats this standard quantum limit constraint when 



C> {l + 2nc)/8. 
Typically the thermal occupation of the cavity, 



(32) 



small, such that the constraint of Eq. ( 32 ) is not exper- 
imentally demanding. Indeed, henceforth, we shall set 
fie — 0. The results accounting for non-zero fie are readily 
obtained by multiplying cavity-dependent contributions 
to the noise spectra by {2fic + 1). 



B. Asymmetric optomechanical couplings and 
mechanical damping 

In the presence of coupling asymmetry, the measured 
collective quadrature is now X^ (rather than A^) as de- 



fined in Eq. (18a), and so it is the noise spectrum of this 



quantity that we seek. Our goal is to understand whether 
coupling and damping asymmetries can be tolerated, or 
even exploited, for the purpose of force sensing beyond 
the usual quantum limits. 

In the presence of asymmetry, the two sets of collective 
quadratures {X+,P-) and {X-,P+) are coupled to one 
another, increasing the complexity of the spectra. In 
general, we write the noise spectral density as the sum of 
a contribution from thermal and quantum fiuctuations of 
the mechanical oscillators, a contribution from the back- 
action of the cavity detector, and an imprecision noise 
contribution, as: 



cmcas r 



w 



SfU+SfM + Sf^iuj]. (33) 



Now in the fully symmetric case considered previously, we 



have Sf [w] = Sx+[uj] as given in Eq. (27), Sf [w] = 



(since the measurement is BAE), and D^^^[a;] is given 



byEq. (31) 



Now the simplest contribution to the spectrum in the 
asymmetric case is the imprecision noise contribution. 



This is given by, c.f. Eq. (31), 



sfn^] = 



1/7 
8(7' 



(34) 



where 7 is now the average mechanical damping rate, 

7 = (7a+7b)/2, (35) 

and we have introduced the rotated cooperativity param- 
eter by 



C = ^^ = C[l + {GdlGf 
7K 



— C sec 



, /arctanp 



(36) 



which follows from Eqs. (20a), (20c) and (29) 



The thermal and back-action contributions to the noise 
spectrum exhibit complicated dependencies on system 
asymmetries. They are given, in a general form, by 
Eqs. (C6) and (CIO), respectively. For weak asymme- 
tries, however, one finds that the measured noise spec- 
trum of X^ still has approximately Lorentzian peaks at 
Lj = ±f2 (not f2), as in the symmetric case. The system 
asymmetries will be quantified through the dimensionless 
parameters p and d, which quantify the degree of coupling 
and damping asymmetry respectively. The parameter p 
is defined in Eq. ([20a), while d is defined by: 



d= {la -lb)/ ha + lb)- 



(37) 



Since the general expressions for these noise contribu- 
tions are complicated, we shall focus on two special cases: 
the noise at the effective mechanical resonance (f2 where 
ri S> 7) and the noise Jar-detuned from the eff'ective me- 
chanical resonance (17 -t- A where $7 ^ A ^ 7). Fur- 
thermore, we shall focus on the limit 7/f7 — > 0. First we 
consider the system in the absence of the compensating 
parametric driving, such that the Hamiltonian is given by 
Eq. ([9]). In Sec. Ill B 3 we shall consider the compensated 
system, with Hamiltonian given by Eq. (25). 



x^ 



1. Noise at mechanical resonance 

First we consider the noise contributions at the ef- 
fective mechanical resonance frequency. The noise con- 
tribution from thermal and quantum fluctuations is a 
complicated function of the coupling asymmetry, as this 
asymmetry modifies the effective mechanical susceptibil- 
ity. However, in the case where the coupling asymmetry 



is zero {p — 0), we simply have 



SfJQ] = -(1/2- 



-(1/2 + n,). (38) 
7& 



That is, the noise contribution from each mechanical os- 
cillator is simply given by the thermal and quantum fluc- 
tuations of the bath to which it is coupled, scaled by its 
damping rate. Another case of interest, for reasons that 
shall become clear later, is that of matched asymmetries 
(p = —d), and in this case we find 



1 



5f [n] = - [(1/2 + fia) (1 + d"") + (1/2 + fib) (1 



(39) 




FIG. 2; (Color online). 





(u - n)h 
Back-action noise, S\ 



[uj], for a 



range of optomechanical coupling and mechanical damping 
asymmetries, centred on oj = +f2. The plotted spectra are 
offset and normalised by the noise spectral density at Q, in 
the ideal, fully symmetric case, c.f. Eq. ( |27[ |. In the fully 
symmetric case (black line), there is no back-action contribu- 
tion to the spectrum. Asymmetries lead to back-action con- 
tributions to the spectra (solid lines). However, by matching 
asymmetries the back-action contribution on resonance may 
be nullified (dashed lines), see Eq. (40 1. The asymmetries 
for each line are (c.f Eq. ([37|): d = 0, Gd/G = (black); 
d = 0.05, Gd/G = 0.025 (blue, solid); d = 0.05, Gd/G = 
-0.025 (blue, dashed); d = 0.2, Gd/G = 0.1 (purple, solid); 
d = 0.2, Gd/G = -0.1 (purple, dashed); d = 1, Gd/G = 0.414 
(red, solid); d = l,Gd/G = -0.414 (red, dashed). The other 
parameters used here (and in subsequent figures) are based 
on the experiments performed by Teufel and co-workers [29] . 
as described in Sec. IV Gl 

Next we consider the noise contribution due to the 
back-action of the cavity on the mechanical oscillators. 
At the resonant peaks, this is 



g.ar^. i(p+^);(i+p;)^. 



X+^ J 7 (1- 



7 (l-d2+p2) 



2^2 



(40) 



That is, the asymmetries lead to back-action heating 
which is proportional to the measurement strength C . 
As shown in Fig.[2| in general, the back-action noise con- 
tribution increases in the presence of either coupling or 
damping asymmetries (as indicated by the solid lines). 
Surprisingly, however, by having appropriately tuned 



damping and coupling asymmetries such that p -|- d = 0, 
one can cancel this extra heating effect at resonance (as 
indicated by the dashed lines) . The measured observable 
is driven by back-action from both the coupling asym- 
metry and damping asymmetry; at resonance these are 
precisely out-of-phase and may coherently cancel. Unfor- 
tunately, though, this cancellation occurs only in a small 
bandwidth about the mechanical resonance, as will be 



described quantitatively in Sec. IV Recall that this back- 
action noise cancellation occurs in the limit 7/5I — > 0. To 
higher-order in 7/ri, and in the case of matched coupling 
and damping asymmetries {p = —d) we have: 



sf p 



7\2 1 



l/7\ 

7 \n) 

x{l + d^f{l + d^-^l + d?^C. (41) 



This result will be useful in determining the ultimate 
limit to force sensing in this system. 



2. Noise away from mechanical resonance 

Next we consider the noise spectrum far from the ef- 
fective resonance frequency (w = ±f2), at a detuning A 
where fi ^ A ^ 7. The thermal contribution in this 
limit is 



s"-^ [n + A] 



7 

4A2 



(1 + d) 



1 



1- 



P 

l+p2 



+ {l-d)[-+m 



l+p' 



(42) 



In the case of no asymmetries (p, d — 0), Eq. (42) reduces 
to a simple result consistent with Eq. (27). 

Next, the back-action contribution to the noise spec- 
trum at w = ±(r2 + A) is 



S'f [f7 + A] = 



_1 

4A2 1 



P 



,C. 



p^ 



(43) 



Clearly, far from the effective mechanical resonance, the 
back-action heating due to the presence of damping 
asymmetry is strongly suppressed compared to the heat- 
ing due to coupling asymmetry. As we will see in Sec. |IV[ 
this will allow for excellent force sensing of a force applied 
to one of the mechanical oscillators which is strongly de- 
tuned from resonance. In the absence of coupling asym- 
metry, the dominant back-action noise contribution is at 
higher-order in 7/ A, 



sfm 



7 / 7 y 

4A2 V2A/ 



d'C. 



(44) 



That is, the back-action heating due to damping asym- 
metry is suppressed by the small factor (7/2A)^. 



3. Noise with compensation 

Now we turn to the noise spectra in the presence of 
the compensation described in Eq. (25 1. The question is 



whether or not such compensation is useful for the pur- 
pose of force sensing, and to answer this question, we 
calculate the measured noise spectra. For weak asymme- 
tries, the noise spectrum remains well-approximated by 
Eq. (27) if we make the replacement Q. ^^ Ct. That is, 



the effective collective mechanical oscillator frequency is 
now VL (rather than fJ), but again we shall quote results 
in two cases: at resonance (f2, assuming f2 3> 7), and far- 
detuned from resonance (f2 -I- A, assuming f2 ^ A ::!> 7). 
At resonance, the thermal noise contribution is 



s'tm 



[Ua -f nf, -I- 1 



7[l-d2(l_p2)] 

+ \nb - riajd 



l-rf2(l_p2) 



(45) 



This reduces to the result of Eq. ( 38 1 in the case p ~ Q 
(where fi — >■ fJ). The back-action contribution is 



st^m^ 



1 



7[l-d2(l-p2 



rC. 



(46) 



Clearly, the back-action heating due to coupling asymme- 
try is now attenuated compared with the uncompensated 
case, c.f. Eq. (|40]). 

Far from the effective mechanical resonance, the ther- 
mal noise contribution, S"*^ [fi + A] , is given by Eq. ( 42 ) 

with p = 0, and the back -act ion noise contribution 



S^ [17 + A] , is given by Eq. ( 44 1 provided that we make 



the replacement C ~^ C. In both cases, the noise contri- 
bution due to coupling asymmetry is suppressed by the 
use of compensation. 



C. Summary 

The noise spectrum of the measured observable may 
be expressed as a sum of contributions from thermal and 
quantum fluctuations, back-action noise and imprecision 
noise, as per Eq. (33). In the fully symmetric case there 



is no back-action contribution, and the noise spectrum 
is given by Eqs. (27) and (31 1. In the asymmetric case. 



the back-action noise contribution is a complicated func- 
tion of the asymmetries. At resonance it is possible to 
cancel this contribution by matching asymmetries (see 
Eq. (40 )), while the contribution far from resonance is rel- 



atively insensitive to damping asymmetry (see Eq. (44|) 



IV. FORCE SENSITIVITY 

A key motivation for the two-mode BAE scheme is the 
possibility of continuously monitoring both quadratures 



of a narrow-band force without any back-action-related 
quantum limit. We now analyze this possibility in detail, 
paying close attention to the role of mechanical dissipa- 
tion, something that was not discussed in Ref. [TUl 



A. Conventional quantum limits 

We consider the standard situation where a narrow- 
band force F{t) is applied to one of the two mechanical 
oscillators in our setup (say the a mechanical oscillator): 



Hr. 



F{t)xa^f{t){a + a^) 



(47) 



where /(i) = /S.XaF{t). In the standard way, we will ex- 
press the total noise in the measured homodyne signal 
as an equivalent added thermal noise on the driven os- 
cillator. The measured force noise spectral density has 
contributions both from the added noise of the measure- 
ment as well as the inherent quantum (zero-point) and 
thermal fluctuations of the measured system, and takes 
the general form: 

S'^'^'^'luj] = hmjau^ (1 + 2na[uj] + 2naddM) , (48) 

where na[w] describes the thermal occupation of the 
driven oscillator and naddii^] is the noise added by the 
measurement. The most straightforward way to moni- 
tor F{t) would be to continuously monitor the position 
Xa of the driven mechanical oscillator. In this case, the 
full quantum limit on continuous position detection [5] 
directly leads to a quantum limit on force detection: 



"■addM > 1/2. 



(49) 



We will refer to this as the full quantum limit on contin- 
uous force detection: the added noise at each frequency 
is at best equal to the zero-point noise. 

Reaching the above-defined conventional quantum 
limit at frequencies far from the mechanical resonance 
requires strong correlations between the position detec- 
tor's back-action and imprecision noises [5], something 
that can be difficult to achieve. In the absence of such 
correlations, one is subject to a more severe constraint, 
the so-called "standard quantum limit" . For frequen- 
cies uj — Ua + A. far from the mechanical resonance (i.e. 
|A| 3> 7a), this standard quantum limit takes the form: 



?^add[Wa + A] > A/7q. 



(50) 



By now calculating the added noise fiadd ['-^] of the two- 
mode BAE scheme which is the focus of this paper, we 
can determine whether one can surpass these conven- 
tional quantum limits. Naively, one might think that 
as one is evading the measurement back-action, there 
should be no quantum limit on the added noise. The 
situation is, however, more complex. The auxiliary os- 
cillator in our scheme (mechanical oscillator b) will have 
its own quantum and thermal fluctuations which con- 
tribute to the output noise and thus the added noise of 
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the measurement. One cannot simply set this to zero by 
making 7^ = 0, as if 7a ^ "ft, the complete BAE nature of 
the scheme is lost (as shown above). Despite these com- 
plications, we find that there are ways to perform force 
sensing using our scheme that yield added noise numbers 
far below the conventional quantum limits. 



B. Results for two- mode BAE 

The two-mode BAE scheme is capable of measuring 
both quadratures of the signal force F{t) as long as it 
is contained in a bandwidth <§; f2 about the mechanical 
resonance frequency oJq ; we thus restrict attention to this 
case. The linearized equations of motion readily yield the 
linear-response relation between the measured observable 
X_|_ and the frequency components of the applied force 
(see App. Id]). For frequencies S satisfying \S\ <C rj, we 
find 

{X+[±iu2+6)]) ^ AXaXF[±ic02+5)]F[±{LJi+S)], (51) 

where 0^2 = ^ and wi — uja in the uncompensated 
case, and XfM is a transfer function (which satisfies 
Xf[— 1^] — Xf[^]*) whose form follows from the equa- 
tions of motion; we discuss its properties further in what 
follows. Note that in the compensated case the rel- 
evant frequencies are shifted such that ^2 = ^ ^nd 
LUi — Lo'^ = uja -\- 0(1 — COS 26*) with 9 as defined in 
Eq. (17), and for this reason results in this section shall 



be expressed in terms of the general frequencies wi and 
UJ2- We shall consider both the cases of a mechanically 
resonant signal force {5 = 0) and a mechanically non- 
resonant signal force (^ = A ^ 7). We stress here that 



the left-hand-side of Eq. ( 51 ) involves the frequency com- 
ponents of Xj^ [lS\ in the rotating frame used to write the 



system Hamiltonian in Eqs. (|9| and ([25|), whereas the 
frequency components of F[Lu\aie in the (non- rotating) 
laboratory frame. 

Eq. (51) tells us that the full spectral information of 



the signal force is reproduced in the dynamics of the col- 
lective mechanical quadrature X^ , the only proviso being 
that the mechanical resonance frequency uji is effectively 
shifted to W2. It follows that the measured force noise 
spectral density is given by: 



S'5?'='^''[cji+(5] = 



\xfV^2 + 5] 



^^'=-^[6^2+ -5], (52) 



where the calculation of the output noise S'^'^^[u!2 



5] 

was discussed in Sec. |III[ recall that these noise spectra 
were explicitly symmetric in frequency. 

Writing the measured force noise spectral density in 
the form of Eq. (48), we can express the added noise of 
our scheme as 



fiaddH = n:^SM + <5dH 



(53) 



where "■addl'^] ^^ ^^^ added noise associated with the in- 
herent thermal and quantum fiuctuations of the auxiliary 
mechanical oscillator, and n'. 



addL 



is the added noise as- 
sociated with the coupling of the mechanical oscillators 
to the cavity. This second contribution may be further 
decomposed into back-action noise, Jiadd['^]i ^^^ impre- 
contributions. These quantities fol- 



cision noise, ^ladd''^ 



low from Eqs. (48) and (52) 



While full expressions for the frequency-dependent 
added noise nadd[w] may be readily derived, for simplic- 
ity we focus on the two main cases of interest: that 
of a mechanically-resonant signal force, and that of a 
far-detuned mechanical signal force. Henceforth in this 
section we consider results in the zero-temperature limit 
{fia, fib = 0), allowing the possibility of force sensing near 
and beyond quantum limits. 



1. Detection of a mechanically resonant force 

The resonant signal force corresponds to having F{t) 
centred in a narrow bandwidth ^ 7^ about the mechan- 
ical resonance frequency ui. We can consider resonant 
force detection both wit hout a nd with the compensation 
scheme described in Sec. Ill B From Eq. ( 51 ) , we see that 



we need to know the transfer function xp at the effective 
resonance frequency a;2 (5 = 0). One finds, in the regime 
^2 > 7, that: 



XF[±i^2] 



T-9r{p,d.), 

7 



(54) 



where gr describes the modification of the susceptibility 
due to asymmetries; its full form as a function of the 
asymmetry parameters p and d is given by Eq. ( D5 ) or 



(D6) in the uncompensated or compensated cases, re- 



spectively. 

In the perfectly symmetric case (i.e. both mechanical 
oscillators have identical optomechanical couplings and 
damping rates), gr — I, implying that X+ responds to 
the applied force analogously to an oscillator driven on 
resonance. In this case, from Eq. (31) we find that the 



added noise of the force sensing scheme on resonance is 



r 1 1 
J^addKJ = ;:: + 



1 

8C 



(55) 



As expected, there is no back-action contribution to 
the added noise, and the added noise is a monotoni- 
cally decreasing function of the measurement strength 
C. Nonetheless, one cannot beat the standard quantum 
limit in this case. Even though the measurement back- 
action goes to zero, the added noise is still limited by 
the zero-point fluctuations of the auxiliary mechanical 
oscillator. 

Now we consider force sensing in the presence of asym- 
metries. First we consider the case of damping asymme- 
try without coupling asymmetry (d 7^ 0,p = 0). In this 



10 



case the added noise due to the auxihary mechanical os- 
cillator is simply 



1 1 



'addl 



21-d 



(56) 



This contribution may be understood in the following 
manner. As the damping rate of the auxiliary mechanical 
oscillator is increased beyond that of the driven oscillator 
(i.e. d from — > —1), its noise spectrum is broadened. 
Therefore most of its noise is outside the force detection 
bandwidth of the driven oscillator, and the added noise 
of the force detection scheme at resonance is reduced. 
If we consider the opposite limit, in which the damping 
rate of the auxiliary mechanical oscillator is reduced be- 
low that of the driven oscillator (i.e. d from — )■ 1), the 
added noise contribution is increased. The force noise 
spectrum due to the auxiliary oscillator is narrowed, be- 
coming sharply peaked about resonance. In either case, 
however, the damping asymmetry alone will lead to a 
back-action noise contribution, as given by Eq. (401, that 



will typically dominate the noise due to the auxiliary os- 
cillator. 

Next we consider the case of coupling asymmetry with- 
out damping asymmetry {p ^ 0,d = 0). Again, coupling 
asymmetry leads to back-action heating and therefore 
reduces force sensitivity. However, compensation, as per 
Eq. (25 1, can be used to almost restore the system's force 
sensing capability. The added noise at the shifted me- 
chanical resonance frequency, w^, is now 



'^addKJ 



1 , / arctannN 1 1 
- sec^ — -^ H ^ ^ - , 

2 V 2 y 8C 2' 



(57) 



with the limit (equal to the full quantum limit) being 
approached in the high-cooperativity, low coupling asym- 
metry regime. Again, the full quantum limit cannot be 
surpassed, however. 

The situation becomes more complex with both damp- 
ing and coupling asymmetries. In general, these asymme- 
tries cause back-action heating and reduce the system's 
force sensitivity. Indeed, using Eq. ([40|), the added noise 



due to back-action, to second-order in p and d, is 



-ba 
'■add 



[LOa]^{p + dfC. 



(58) 



The added noise now has terms which increase with the 
measurement strength C. 

However, from Eq. ( [40| , we know that it is possible to 
cancel back-action on resonance by matching asymme- 
tries (i.e. setting p = —d). This is promising for force 
sensing, but we must also check the added noise due to 
the auxiliary mechanical oscillator. This contribution is 



nlTA^a] 



1 1 



21 



-{l + d + d^){l + 



VTTd^ 



(59) 



Clearly, the added noise contribution of the auxiliary me- 
chanical oscillator goes to zero when d — >■ 1, and com- 



bined with Eq. (401, this suggests the added noise of the 



scheme can go to zero with matched asymmetries. This 
suggestion is based on calculations performed in the limit 
7/n — >■ 0. We can check this possibility by calculating 
the added noise to higher-order in 7/ri. The added noise 
contribution due to back-action with matched asymme- 
tries {p = —d) is 



^addl^a 



1 - Vl + rf2 -d{l-d- Vl-(i2) 



(1 + ^2)3 1 ,^ 



1 + d 8 



(g)- 



C. 



(60) 



asymmetry case {d 
due to the cavity is 



From Eqs. (41 ) and ( |60[ ) we find, in the extreme damping 
1 andp 



'addl 



rC 



-1), that the added noise 



(61) 



21^2 128(^/2 -1)C' 



While there is a back-action contribution here (first 
term), it is proportional to (7/^7)^ and hence small. One 
readily finds an optimal measurement strength. 



Co^ 



V3 n 



17 



at which point the total added noise is 



"add[WaJ 



Co 



7 



V3 



1 



(62) 



(63) 



The added noise remains close to zero provided that 
57 ^ 7, and so the full quantum limit can be beaten 
for resonant force sensing via matched asymmetries. 

Unfortunately, this result applies over only a very small 
bandwidth about the mechanical resonance frequency. 
The added noise may be evaluated as a function of fre- 
quency, and we are interested in the bandwidth B over 
which the added noise remains below the full quantum 
limit. It may be shown that, in the case of extreme damp- 
ing asymmetry (d — >■ 1), this bandwidth is 



B 



1_ 7 
2V2VC 



7 



(64) 



with the second scaling attained for an optimal measure- 
ment strength (cooperativity), see Eq. (62). That is, the 
useful bandwidth becomes very small. 

Plots of the added noise, as a function of the damping 
rate of the auxiliary oscillator (with the damping rate of 
the driven oscillator held constant), are shown in Fig. Is] 
The added noise is shown for three coupling asymmetries, 
and in each case the added noise at resonance may be 
reduced to the level of the auxiliary oscillator fluctuations 
by matching coupling and damping asymmetries. 



2. Detection of a mechanically non-resonant force 

We now turn to the case where the signal force F{t) 
is detuned from the mechanical resonance, and again we 
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FIG. 3: (Color online). Added noise for resonant {u — uja) 
and non-resonant {uj — iUa + A) force sensing, nadd[i.tj], as a 
function of the damping rate of the auxiliary mechanical os- 
cillator, 7i,/7a, while 7^ is held fixed. The dashed black line 
is a lower bound on the contribution due to auxiliary oscilla- 
tor thermal fluctuations; see Eqs. ( 59 1 and ( 68a I , for example. 



The blue, purple and red curves correspond to the added noise 
for resonant force sensing at Gd/G = 0,-0.162 and —0.414, 
respectively. The back-action contribution to the added noise 
can be made to go to zero on resonance by matching asymme- 
tries, irrespective of the size of the damping asymmetry. The 
cyan line is the added noise for off-resonant force sensing, 
with no coupling asymmetry and the measurement strength 
optimised. 

will consider both the original and compensated systems. 
Specifically, we consider the case where F{t) is contained 
in a bandwidth B centred on a frequency uji + A, where 
|A| » 7a and B < |A|. From Eq. (51 1, we see that we 
need to know the transfer function at ^2 + A, 



1 



Xf[±(w2 + A)] = ^9niP,d), 



(65) 



where (?„ describes the modification of the transfer func- 
tion due to asymmetries; its full form as a function of 



the asymmetry parameters p and d is given by Eq. (D7) 



in the original case and by Eq. ( D8 ) in the compensated 
case. As compared with the case of a mechanically res- 
onant signal force, see Eq. (54), the symmetric "gain" is 



now 1/(2 A) rather than I/7. Also note that force detec- 
tion, in this case, is not limited by a resonant detection 
bandwidth of order 7. 

In the perfectly symmetric case (i.e. both mechanical 
oscillators have identical optomechanical couplings and 
damping rates), g„ ~ 1. Then we find that the added 
noise is 



nadd i^a 



Al = 



1 



t) 



1 

8C 



1 
2' 



(66) 



with the limit being approached in the high-cooperativity 
regime. Clearly, the scheme surpasses the standard quan- 
tum limit on non-resonant force detection (c.f. Eq. ([50|). 
It, however, is only equal to the full quantum limit (c.f. 
Eq. (49)) that comes from allowing detector noise corre- 



lations. Again, the residual added noise of 1/2 a quantum 



in Eq. ( 66 ) is due to the zero-point noise of the auxiliary 



mechanical oscillator. 

We can easily determine the effects of asymmetries on 
our force sensing scheme. From Eq. (43), to second-order 
in the asymmetry parameters and in the limit 7/A — >■ 0, 
we find 



'^addl 



A] 



p'C. 



(67) 



As already noted, in this highly-detuned regime, it is only 
the optomechanical coupling asymmetry that yields an 
appreciable deviation from being quantum-limited. This 
suggests that one could improve on the symmetric added 
noise result of Eq. ( 66 ) by exploiting damping asymmetry 



to reduce the added noise contribution from the auxiliary 
oscillator. 

Indeed, these added noise contributions due to the aux- 
iliary mechanical oscillator are given by 



«S3K + A] 



1 l-d 



21-l-d 



i+p{p+ yrTp2) 



'add 



K + A] = 



1 1 

2Y 



d 



: sec 



/arctanp 
V 2 



(68a) 
(68b) 



in the original and compensated cases, respectively. In 
both cases, they go to zero in the limit d — >■ 1. Combined 
with the insensitivity of the back-action noise to damping 
asymmetry as expressed in Eq. (67), this would appear 



to enable force sensing beyond the full quantum limit, 
provided that there is no coupling asymmetry. In fact, we 
might expect that this force sensitivity can be achieved 
in the presence of coupling asymmetry if compensation 
is used. 

These arguments can be made precise. The following 
results describe the case with compensation included, and 
so reduce to the case of no compensation when there is 
no coupling asymmetry to begin with {p = such that 
o;^ — > oJq). The added noise associated with the cavity is 



'add l^a 



7^ d^ A /arctanp 

—tt; sec 

4A2 1 + d V 2 



C 



+ ■ 



4A2 1 1 
^l + d8C' 



This added noise is minimised by a cooperativity. 



Co = 



4A2 



1 



cos 



/arctan p 



7' 2y2|d| ""^ V 2 
at which point the added noise due to the cavity is 

1 Ml 



(69) 



(70) 



^aSdK + A]! 



Co 



^|l + d| V 



o / arctan n\ 
sec^ I ^— ^ 1 



(71) 



Taking the limit d — > 1, the added noise due to the aux- 



iliary oscillator vanishes as per Eq. (68b), and the total 
added noise becomes 



naddK + A]|^ = 



1 / arctan bN 1 

— - sec^ — -^ -^ — ^ 

2V2 V 2 y 2^/2 



, (72) 
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with the result tending to the limit as the coupling asym- 
metry, p —>■ 0. Now, one can not only beat the stan- 
dard quantum limit for non-resonant force detection (c.f. 
Eq. ( |50[ )) by a greater amount, but one can also surpass 
the full quantum limit on force detection (c.f. Eq. (49 1). 
Without compensation, the added noise in Eq. ( 72 1 grows 
more rapidly as \p\ increases. 



C. Summary 

In summary, for sensing of a mechanically-resonant 
force, force sensing at the full quantum limit (being equal 
to the standard quantum limit, on resonance) is possible 
in the symmetric case. Even with coupling asymmetry, 
the possibility of force sensing near this quantum limit is 
retained via compensation. Further, it is possible to sur- 
pass this full quantum limit via matched asymmetries, 
though only in a very narrow bandwidth. 

For sensing of a mechanically non-resonant force we 
can surpass the standard quantum limit in the symmetric 
case. Furthermore, one may exploit damping asymmetry 
to surpass the full quantum limit. In both cases, these 
results can be retained in the presence of coupling asym- 
metry using compensation. Additionally, force sensing 
in this case is not limited to a resonant detection band- 
width. 



the single-phonon levels have both been demonstrated, 
though this is not the case for mechanical continuous 
variables. The entanglement of collective spin operators 
of atomic ensembles [3F, and of the motional states of 
trapped ions [36^ . have both been achieved, however. 

There has been a great deal of discussion of the pos- 
sibility of entangling macroscopic mechanical degrees of 
freedom with other degrees of freedom (see [35] and ref- 
erences therein), including other mechanical degrees of 
freedom. Particular attention has been paid to the possi- 
bility of all-mechanical entanglement with the interaction 
being mediated by an electromagnetic field. Proposals 
studied include placing mechanical oscillators in a ring 
cavity or interferometer, both without [39ti4T] and with 
|42j measurement and feedback, placing dielectric mem- 
branes in a cavity [43^ , and using remote optomechanical 
systems, both with [44ti46] and without [471 148] protocols 
dependent on optical measurements. 

More recently, Schmidt and co-workers have proposed 
an all-mechanical entanglement generation scheme based 
on the detuned, modulated (and therefore, two-tone) 
driving of a coupled electromagnetic cavity |49) . Com- 
pared with their approach, our scheme has the consider- 
able advantages that it does not require highly-detuned 
driving of the cavity, and is robust to the initial thermal 
populations of the mechanical oscillators. Further, the 
possibility of entanglement verification is built into our 
proposal. 



CONDITIONAL VARIANCES AND 

ENTANGLEMENT 



B. Conditional Dynamics in Adiabatic Limit 



In Sec. |III| we calculated the noise spectra of measured 
and perturbed observables in our system, with a view to 



analysing its force sensing capabilities in Sec. IV These 
are all unconditional quantities; that is, the description 
of the state of the system is not updated based on the 
measurements made upon it. However, since we are con- 
tinuously monitoring the system, we can also discuss its 
conditional dynamics, arising from updating the best es- 
timate of the state of the system based on the measure- 
ment record. In particular, a continuous measurement 
conditionally projects the system into an eigenstate of 
the measured observable. This shall enable, in this sys- 
tem, conditional mechanical two-mode squeezing and the 
entanglement of the two mechanical oscillators. 



A. All-Mechanical Entanglement 

The preparation and verification of macroscopic, all- 
mechanical entanglement is a fundamental goal of the 
study of mechanical systems in the quantum regime [7] . 
Such a state is also a physical approximation of an 
Einstein-Podolsky- Rosen channel [3l, a key ingredient 
in quantum information processing protocols with con- 
tinuous variables [32. Electromechanical entanglement 
[55] and the entanglement of phonons in bulk [33], at 



The conditional dynamics of the system may be de- 
scribed in the standard manner using a stochastic master 
equation [SD|. The evolution of the joint density opera- 
tor, (J, of the two mechanical modes and one electro- 
magnetic mode, under homodyne detection of the cavity 
mode quadrature ce*'^-|-c^e~*"^ with a quantum efficiency 
r], is given by 

da = —i\H, cr] dt -f- K{nc + 1)2? [c] a dt 

+KnJ) [c'f] a dt + y/^M[ce"'']a dW 
+la{na + l)'D[a\a dt + "fana'D[a^a dt 
+jb{nb + l)V[b]cr dt + jbnb'D[P]a dt, (73) 

where H denote s th e ori gina l or compensated Hamil- 
tonian of Eq. (19) or (25), respectively, T>[A]a = 



AcrA^ — ^A^Aa— ^aA^A is the dissipative superoperator, 
A^[A]cr = Aa + a A^ — Tr [Aa + aA^a is the measurement 
superoperator, and dW is the Wiener increment. 

We consider the "good measurement" limit, in which 
the cavity damping rate exc eeds the rates in the Hamil- 
tonians of Eq. (19) or (25); that is, k > fl,pri.,G (f^ 



rather than fi in the compensated case). In this limit, the 
cavity responds rapidly to the dynamics of the coupled 
mechanical oscillators, and we may adiabatically elim- 
inate the cavity mode [5T]. Setting (f> — tt/2 for con- 
venience, the cavity mode lowering operator is given by 
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c = —i\/2GXj^j K. Note that it is crucial that the adi- 
abatic eUmination is performed in terms of the rotated 



observables; otherwise, the dissipative terms in Eq. ( 73 ) 



will involve linear combinations of mechanical annihila- 
tion operators, which greatly complicates the description. 
Now the evolution of the (reduced) joint density oper- 
ator of the two mechanical oscillators having traced out 
the cavity mode, p = Trcav [c], is given by 



In general, the elements of S will describe covariances 
of the rotated observables, though we drop the tilde no- 
tation in the subscripts of the off-diagonal covariances 
in Eq. ( 78 ) for convenience. In the symmetric case, the 



tilde notation on the subscripts of the collective quadra- 
ture variances in Eq. (78 1 shall be dropped. 



dp 



The equations for the best estimates of the collec- 
tive (sum and difference) mechanical quadrature observ- 
ables take the form of an Ornstein-Uhlenbeck process (c.f 
Eq. (Ill), 



-i\H! ,p\ dt 

-^[X+, [X+,p]] dt+^M[X+]p dW 

-la{na + l)'D[a\p dt + -^anaV[a}]p dt 

-7b(nfc + l)V[b]p dt + -fbfn,V[P]p dt, (74) where the system matrix M is given by Eq. (|B9|) or ( |B12| ), 



dt 



F = M-y + Q-e(t), 



(79) 



where %' now denotes either the Hamiltonian of Eq. ( 19 ) 
or (25), excluding terms involving the electromagnetic 
mode operators. We have also introduced the collective 
quadrature measurement rate as 

r = 7C' = 2&/k, 



(75) 



where G is the modified coupling constant introduced in 
Eq. ( 20c ) . The measurement record increment associated 



with Eq. ( 74 ) is given by 



dr = {X+)dt + dW/y 



(76) 



Integrating Eq. ( |74[) a ccording to the measurement record 
increment of Eq. ( |76[ ) allows us to continually update our 
best estimate of the state of the two mechanical oscilla- 
tors. 



Best Estimates of Quadratures and Conditional 
Variances 



Now according to Eq. ( 74 1 , with the quadratic Hamil- 
tonian obtained from Eq. (19 1 or (251, and linear damp- 
ing, decoherence and measurement, the steady-state re- 
duced density operator p will be Gaussian. Therefore, 
the conditional state of the system will be fully described 
by knowledge of its first and second moments. One can 
immediately write down the equations of motion for the 
(scalar) best estimates of the collective quadrature ob- 
servables and for the conditional variances corresponding 
to the stochastic master equation in Eq. (74 1. The vector 
of best estimates is formed as 



y- x+,p_,x_,p+ 



(77) 



For later convenience, we introduce the following nota- 
tion for the ten independent elements of the two-mode, 
symmetrically-ordered covariance matrix, S: 



(78) 



r^x. 


S+- 


^xx 


S++ 


S+- 


Vp_ 


E— 


T,pp 


Sxx 


E__ 


Vx 


S-+ 


Ls++ 


Yipp 


S-+ 


^P. 



for the Hamiltonian of Eq. (19) or (251, respectively, the 



measurement noise weighting vector is 



Q 



^;f,,S+_,E 



XX, 



(80) 



and ^{t) = dW/dt is a white noise process describing the 
noise in the measurement. 

The equation for the best estimates of the covariances 
of the collective mechanical quadrature observables is 



S = MS + SM^ + L - SK^KS, 



(81) 



where S is the symmetrically-ordered two-mode covari- 
ance matrix in the ordered basis defined by E q. (|21| , 
{X+ P^,X^,P+), and again M is given by Eq. (B9[) or 
(B12). The matrices K and L are defined by 



(K)i 



0, 



n 








■"eq 





■"cq 





"oq 



(82a) 
(82b) 



where we have introduced the notation n' = nth + 1/2 + 
dfifi and n'^ — n^ + d(nth + 1/2), in terms of the dimcn- 



sionless damping asymmetry introduced in Eq. (37) and 
the effective mechanical thermal occupations. 



nth 



n-d 



{ua + rib) 
{fia - fib) 



(83a) 
(83b) 



Note that an equation of the form of Eq. (81 1 has previ- 



ously been obtained in a study of the optimal control of 
two-mode entanglement (parametrically-coupled optical 
modes) under feedback [37]. 



Due to the rules of the Ito calculus, Eq. (81 ) is a deter- 



ministic system and we may solve for the steady-state co- 
variance matrix. The steady-state equation correspond- 
ing to Eq. ( [sTj ) is a continuous-time algebraic Riccati 
equation (CARE). This is a system of nonlinear algebraic 
equations, though one for which numerical methods are 
well-developed. 
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D. Steady-State Conditional Variances and 
Entanglement 

The two-mode, symmetrically-ordered covariance ma- 
trix S enables a full characterization of the entanglement 
and purity of the system |52| . Indeed, the entanglement 
may be directly quantified using the logarithmic nega- 
tivity [53^ . From an experimental perspective, however, 
the estimation of all elements of the covariance matrix is 
challenging. A simpler method for determining whether 
or not the mechanical oscillators are entangled is pro- 
vided by Duan's inseparability criterion [28]. Note that 
this is a sufficient, but not necessary, condition for the 
inseparability of a bipartite quantum state. Specifically, 
we consider a generalized version of Duan's criterion |54j . 
that here takes the form 



^x+ + ^p_ < cos 29 






1 



(84) 



In the case of symmetric optomechanical coupling (G^ — 
0), the criterion of Eq. ( [84| reduces to the standard ver- 
sion of Duan's criterion. 



Vx^ + Vp^ < 1. 



(85) 



Crucially, both quantities on the left-hand-side of 
Eq. (84) may be obtained in a straightforward manner 



from the measurement record. The quantity X^ is mea- 
sured directly, while in the regime i7 ^ 7, P-_ is dy- 
namically coupled to it such that it is also effectively 
measured. 



E. The Symmetric Case 

1. Conditional Variances 

In the ideal, symmetric case, an analytic form for 
the steady-state covariance matrix can be found (see 
App.[E|. Of particular interest is the variance of the mea- 
sured collective quadrature. Taking the limit j/fl — ^ 
while the cooperativity C is fixed leads to 



V: 



X4 



V877G(nth + 1/2) + 1-1 V^th + 1/2 1 



ArjC 



277 



C 

(86) 

with the latter result following for C ^ l/[8?7(nth + l/2)]. 
This is a very good approximation for most experimen- 
tally relevant scenarios, including that we shall consider 
in more detail below. Alternatively, taking the limit 
C — 7- 00 while 7/17 is fixed, we find the form 



V477C(nth + 1/2) + 1-1 V^th + 1/2 1 



4r]C 



2v^ 



C 

(87) 

The scaling with C is now the same as in the case 

57 = (c.f. Eq. (E3a)), which itself scales as in the 



case for single-mode back-action-evading measurement 
[6]. Clearly, for a sufficiently strong measurement (large 
cooperativity), the conditional variance of the measured 
collective quadrature observable, X+, is squeezed below 
the vacuum level (1/2) and tends asymptotically to zero. 
To assess whether we are truly generating a two-mode 
squeezed state and entanglement, we also need to check 
that the P_ quadrature is squeezed by the measurement 
(c.f. Eq. (85)). It can be shown, see App. [e| that 



Vp_ = Vx. 



1 



[idT^^C^nto 



l6riC 

+ (1 + i-qCntot + Vl + 8r/Cntot) ' 



to second-order in j/fl, where we have introduced the 
notation for the total thermal and quantum fluctuations. 



ntot = nth + 1/2. 



(89) 



Note that Eq. ( p8| does not imply that the difference be- 
tween Vp_ and Vx^ diverges as C — >■ 00. In this limit 
higher-order contributions are important, and the pertur- 
bative result of Eq. ( [SS] ) is not valid. However, provided 
the ratio j/il is small, from Eq. (88) we do have 



Vp 



Vx, 



(90) 



That is, in the limit that the collective oscillator fre- 
quency greatly exceeds the mechanical damping rate, the 
collective quadrature P_ that is dynamically coupled to 
the directly measured observable X+ is also effectively 
measured and therefore conditionally squeezed. Accord- 
ingly, the measurement conditionally generates a two- 
mode squeezed state. 
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200 300 
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FIG. 4: (Color online). The Duan quantity, Vx, + Vp_, as 
a function of the cooperativity C, for effective temperatures 
corresponding to fith ~ 0,1,25, and in the regime f2 3> 7. 
The bound on the Duan criterion, as in Eq. (851, is shown as 
the dashed black line. The quantity Vx, + Vp_ is seen to be 
below this line for sufficiently large cooperativities, indicating 
that the mechanical oscillators are entangled. This remains 
the case even for a relatively large thermal occupation. The 
required cooperativities are experimentally accessible. 

That both A"+ and P_ are effectively measured (and 
therefore conditionally squeezed) in the regime fi 3> 7 
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can be seen in the following way. Defining primed oper- 
ators by the transformation Zj_ = ^-iHot ^_^^+iHot t^[^Ij 
Ho = ^{X+X- +P+P-), we can express the observables 



i.oor 



_i_ and P_ as 








\x+] 




cos Qt sin V,t 




P_ 




— sin fit cos fit 



P' 



(91) 



While the observables X^ and P_ oscillate at ±51, the 
observables X^ and Pi are constants of the motion with 



respect to the closed system dynamics. From Eq. (91) 



it is clear that by continuously monitoring X_|_, we are 
continuously monitoring both X'^ and PL- Since these 
observables also determine P_, we are effectively mea- 
suring P_ as well (provided that $7^7). 

Recalling the "subspaces" (or subsystems) introduced 
in Sec. |llj we now know the conditional variances of the 
observables in the "measured" subsystem. For the ob- 
servables in the "perturbed" subsystem, the variances in 
the hmit j/fl — > follow from Eqs. (E9a) and (E9b) as 



V- 



X- 



•Vp^ 



nth + 1/2 + C/2. 



(92) 



Effectively, we are measuring both X^ and P_ in this 
limit, and so the conjugate observables of both are 
equally perturbed. 



2. Entanglement 

Next we explicitly consider the entanglement of the 
mechanical oscillators. From Eqs. (86) and (90), in the 



regime fi 3> 7 and the strong measurement limit, we 
find that Vx^ + Vp_ < 1, implying that our mechanical 
oscillators are conditionally entangled. Substituting the 
asymptotic form of Eq. (86) into Duan's criterion leads 



to the sufficient condition on the measurement strength 
for the generation of mechanical entanglement by mea- 
surement. 



C> 



2(nth + 1/2) 



V 



(93) 



Experimentally, this is not an overly demanding condi- 
tion. It is interesting to note, however, that this condi- 
tion does not exhibit the total insensitivity to tempera- 
ture that has been found in a scheme to generate an EPR 
channel in a cascaded atomic and nanomechanical system 
[13j . The primary distinction is that here we are describ- 
ing steady-state entanglement, rather than instantaneous 
entanglement arising from a strong single-shot feedback 
operation. 

The Duan quantity of Eq. ( 85 ) is plotted in Fig. l4J as a 
function of cooperativity, for a range of thermal occupa- 
tions and for parameters corresponding to the experiment 
of Teufel and co-workers [29] . It is seen that the mechan- 
ical oscillators will be entangled in an achievable param- 
eter regime, even for mechanical oscillators initially in a 
thermal state far from the ground state. 




K^^^ 0.25 



0.00 0.15 0.30 0.45 0.60 0.75 
Coupling asymmetry, Gd/G 

FIG. 5: (Color online). The variances Vp_ (dashed lines), 
and the generalized Duan quantities Vx +Vp_ (solid lines) as 
a function of the optomechanical coupling asymmetry Gd/G, 
for the original (Eq. ( 19 1) and compensated (Eq. ( 25 1) Hamil- 



tonians. We assumed 7a = 7b, though even for large damping 
asymmetries, the modification to these lines is small. The 
bound on the generalized Duan inequality, as in Eq. ( |84[ ), is 
shown as the curved dashed black line. This plot was pre- 
pared with the effective occupation set to nth = and a co- 
operativity of C = 100. At such low temperature, the scheme 
is highly robust against coupling asymmetry. Even without 
compensation, the mechanical oscillators remain entangled for 
a coupling asymmetry of up to ~ 60%. With compensation, 
the scheme is unaffected by coupling asymmetry. 

Further, one can explicitly evaluate the entanglement 
measure known as the logarithmic negativity [5 3) . Using 
the result of Eq. ( 86 ) we find 



E^^\ {log2 [riC/inti. + 1/2)] - 1} , (94) 



in the regime specified by the criterion of Eq. (93), and 



zero otherwise. The onset of entanglement, as a function 
of cooperativity, according to this measure is found to 
be consistent with that predicted by the application of 
Duan's criterion (c.f. Eq. (93)). 



F. The Asymmetric Case 

If the damping rates or the optomechanical coupling 
rates of the mechanical oscillators are different, an an- 
alytical solution for the steady-state covariance matrix 
is no longer possible and we must resort to numerical 
methods. The one exception to this is the case where we 
have coupling asymmetry alone and compensation as per 
Eq. (|25|. In this case, the asymptotic results of Eqs. (86) 



and ( |87[ ) remain valid, implying entanglement can still be 
achieved for strong measurements; see App. [E] for more 
details. 

Figs. |5] and [6] show the steady-state conditional vari- 
ance Vp_ and the generalized Duan quantity V^ -|- Vp_ 
as a function of the coupling asymmetry Gd/G. These 
are plotted assuming 7^ = 7^, though the entanglement 
generated is only very weakly dependent on damping 
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asymmetry, even when this asymmetry is large. Curves 
are shown for both the original and compensated cases 
in both figures, for the effective occupation fith = in 
Fig. [5] and for nth = 5 in Fig. [6j Note that in the original 
(uncompensated) case, the variance of the directly mea- 
sured observable V-^ (the difference between the solid 
and dashed lines) is independent of the coupling asym- 
metry, while the variance of the dynamically coupled ob- 
servable Vp is less effectively squeezed as the coupling 
asymmetry is increased. 

At fith = 0, the entanglement generation is seen to be 
highly robust against coupling asymmetry, even without 
compensation. The mechanical oscillators remain entan- 
gled up to a coupling asymmetry of ~ 60%. At fith = 5, 
the robustness against coupling asymmetry is reduced. 
In this case, the oscillators only remain entangled for a 
coupling asymmetry up to '^ 17%. It should, however, be 
possible to independently cool the mechanical oscillators 
via auxiliary cavities such that they are initially at low 
thermal occupations [53]. Further, with compensation 
included, the entanglement generation is totally insensi- 
tive to coupling asymmetry. Furthermore, the coupling 
asymmetry for which entanglement generation is possi- 
ble, beyond a threshold value, is not strongly dependent 
on the cooperativity. 
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Coupling asymmetry, Gd/G 
FIG. 6: (Color online). The variances plotted in Fig. [s] but 
now with an effective thermal occupation fith = 5. At this 
higher effective temperature, the robustness of the scheme 
(without compensation) is reduced. The mechanical oscilla- 
tors remain entangled for coupling asymmetries up to ~ 17%. 
With compensation, the scheme is unaffected by coupling 
asymmetry. 



G. Experimental Parameters 

The experimental parameters used to prepare Figs. [2]- 
[7] correspond to a micromechanical membrane in a su- 
perconducting microwave cavity, a design due to Teufel 
and co-workers [IS]. One could conceivably form the 
microwave cavity from two micromechanical membranes, 
forming the three-mode optomechanical system that we 
have studied. However, the two-mode BAE scheme we 
have described could be performed in a number of cavity 



optomechanics implementations. 

For the purpose of numerical calculations, the cav- 
ity resonance frequency is taken to be 8 GHz and the 
mechanical resonance frequencies are centred around 
10 MHz. The cavity damping rate is K/27r = 200 kHz, 
such that the system is operating well into the resolved- 
sideband regime. The mechanical damping rate is 
7/27r — 100 Hz. The zero-point fluctuations of the mem- 
brane are given by Aa; — 4 fm, such that the nominal 
single-photon optomechanical coupling rate in Eq. ^ is 
gi/27T — 200 Hz. Further, suppose the cavity is driven 
such that its steady-state population at the driven side- 
bands is 10^ photons. This corresponds to a (nominal) 
effective optomechanical coupling rate in Eq. Q that is 
G/27T — 200 kHz. These parameters correspond to a co- 
operativity of C = 4 X 10'^. This is a much larger coopera- 
tivity than is actually required for the scheme presented. 
In fact, such an optomechanical coupling rate would in- 
validate the adiabatic elimination performed in order to 
obtain Eq. (FTil). 



For the purpose of calculations here, we take a nominal 
effective coupling G/2tt — 70.7 kHz, corresponding to a 
cooperativity of C = 500. Throughout this work, a quan- 
tum measurement efficiency oi rj = 1 has been assumed. 
Note that, from the perspective of entanglement genera- 
tion, one can compensate for the reduced efficiency via a 



larger cooperativity (c.f. Eq. (93)). These parameters are 
within the limits of validity of our analysis: recall that 
we assumed resolved-sideband operation (coa jUj^^k) in 
order to avoid spurious back-action heating of the mea- 
sured observable, and the "good measurement" regime 
(k > il.,pCt, G) such that it was possible to perform the 
adiabatic elimination. 



H. Summary 

The conditional dynamics of the system under homo- 
dyne detection of the field output from the cavity are 
described by Eq. (73 1. In the adiabatic limit, in which 



the cavity responds rapidly to the mechanical motion, we 
have the simplified description of Eq. ( 74 ) . This equa- 



tion allows us to calculate the steady-state conditional 
variances of the collective mechanical quadratures (see 
Eq. (81)), and then assess the entanglement of the me- 
chanical oscillators using Duan's criterion (see Eq. (84)). 
It is found that the mechanical oscillators will be entan- 
gled for reasonable experimental parameters (see Fig.|4|, 
even in the presence of considerable asymmetries (see 
Figs. [5] and [6]). 



VI. UNCONDITIONAL VARIANCES AND 
FEEDBACK 

The variances and covariances calculated in the previ- 
ous section as solutions of Eq. (81 ) are conditional quan- 
tities. However, the measured output noise spectrum of 
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Eq. (31 1 is an unconditional quantity. This means it con- 



sists of both the fluctuations described by the best esti- 
mate of the conditional variance (V^ ) and the fluctua- 
tions in the best estimate of the observable itself ((X^)): 






X4 



Vx^+{XI) 



measured noise spectral density, or equivalently, from the 
measurement record. The unconditional variance of X^ 
is obtained directly, while the unconditional variance of 
P_ follows from looking at the quadratures of the mea- 
surement signal itself, as described around Eq. (91). 



(95) 



The fluctuations in the best estimates of X+ and P_ 



may be obtained from the filter of Eq. ( 79 ) . Written out 



We can define the total (unconditional) variances of 
the other collective quadratures in a similar fashion. 
From Eq. (95), the fact that the conditional variance is 



explicitly in terms of the measurement record of Eq. ( 76 ) , 
the filter is 



squeezed below the vacuum level does not necessarily im- 
ply that the unconditional variance will be also. Here we 
calculate the total variances, and describe how one can 
use feedback to reduce the unconditional variances to the 
conditional variances. 



A. Unconditional Variances 

Fluctuations in the best estimates of the collective 
quadrature observables are described by Eq. (79). The 



— V = M-V- 
dt 



^ [i{t) - x+\ g. 



(99) 



where M, Q and V are as specified following Eq. (79), 



and I{t) — dr/dt is called the measurement current. Of 
course, the filter of Eq. (99 ) is itself dependent on knowl- 



edge of the system's covariance matrix. We assume the 
system covariances have reached their steady-state val- 
ues: that is, one should use the calculated values that 



follow from the steady-state solution of Eq. (81) 



steady-state fluctuations are then given by the solution 
of the Lyapunov equation. 



In the fully symmetric case, the measured subsystem 
decouples from the perturbed subsystem. Then the fllter 
of Eq. (|99| can be recast, in the limit S^ — ?► 0, as 



MS + SM^ = -QQ^ 



(96) 



X+{t) 
P. it) 



UriT Vx+e 



-7* 



cosQ.t 
— sin fii 



*/(i), (100) 



where S is a matrix of second moments of best estimates 
of collective quadrature observables, described in the 
ordered basis defined by Eq. (|77|), (X+,P_,X_,P+)^. 
The matrix M is given by Eq. (|13p , (|B9|) or (|B12| , in measurement record of Eq. (|100[) 



where the asterisk denotes convolution and we have in- 
troduced the notation 7 = ■y + 4j]rVx^- The fluctuations 
in the best estimates may be obtained from the filtered 



the symmetric, uncompensated asymmetric or compen- 
sated asymmetric cases, respectively, while Q is given by 
Eq. (15) in the symmetric case and by Eq. (B3) in both 



asymmetric cases. 

In the fully symmetric case, and in the relevant regime 
ri 3> 7, we find the total unconditional variances to be 



C. Feedback 



1. Unconditional Variances 



y^°t 



vl""- = V- 



' X 



x. 



(1 -f 2i^CVx, 



(97) 



Substituting the first asymptotic form of Eq. ( 86 ) into 
Eq. Wi\ we find that 



vk°^ 



t4°* 



nth + 1/2. 



(98) 



This is precisely as expected; a back-action-evading mea- 
surement should neither heat nor cool the mechanical 
oscillators in the unconditional sense. The unconditional 
variances in the asymmetric case, both without and with 
compensation, are given in App. [FT| Note, however, that 
Eq. (97) remains valid (in terms of rotated observables) 
for the compensated case with no damping asymmetry. 



B. Measuring Conditional Variances 

The conditional variances may be obtained from mea- 
surements of the unconditional variances and the fluctu- 
ations in the best estimates, as per Eq. (95). The un- 



As noted earlier, feedback may be employed to reduce 
the fluctuations in the best estimates of our observables, 
and so reduce the unconditional variances to conditional 
variances. In general this may be achieved by adding 
the damping terms — (a7/2)X+ and — (a7/2)P_ to the 
equations for X^ and P^ (respectively) to the filter of 
This requires the application of (asymmetric) 
jk forces to the two mechanical oscillators, given 



Eq. (99) 



feedbac 
by 



Pa(b){t) 



07 
71 



(cos 6 ^ sin 6) X^ sin uj^t 



+ ^ (cos6'±sin6')P_cosa;™i, (101) 
v2 



in the laboratory frame. From the form of Eq. ( 99 ) , we 
expect that X+ and P- will themselves oscillate at ±J7, 



and therefore the applied feedback force of Eq. ( 101 ) will 



conditional variances may be obtained directly from the 



in fact have significant weight at the two mechanical res- 
onance frequencies. In most cases (except the case where 
there are significant coupling and damping asymmetries) 
it suffices to apply the feedback component proportional 
to the estimate X+ alone. 
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In order to calculate the variances under feedback we 
can use an equivalent classical description of the dynam- 
ics of our quantum system [SS] , coupled to the filter giv- 
ing the best estimates of the collective quadrature ob- 
servables (and upon which the feedback will be based). 
The classical quantity representing the collective mechan- 
ical quadrature X^ is denoted x+ (and similarly for the 
other collective quadratures). The total unconditional 
variance of the measured observable, under feedback, is 
theny?' ={xl). 

Now we can write the measurement record increment 
in terms of such a classical representation as dr — x+dt + 
dW/ \/AT]r. We can also write the measurement record 
increment in terms of the best estimate of the measured 
observable as dr = X^dt + dW/^/ArjT. Comparing these 
expressions leads to 



dW = y^ {x+ - X+) dt + dW. 



(102) 



Putting Eq. ( |102[ ) together with the system dynamics 
and filter from Eqs. (12) and (99), respectively, both the 



conditional and unconditional dynamics can be described 
by the Ornstein-Uhlenbeck process, 



dt 



r = -s-y + T- 



dW 

~dF' 



where the system matrix S and the noise weighting ma- 
trix T are defined in Eqs. (|F2]) and (|F5]) of App. [f] 
The state vector Y and noise increment vector dW in 



Eq. ( 103 ) are formed as 



Y = 



V 



dW ^ 



dW 



dw 



(104) 



where V is the vector of best estimates of collective 
quadrature observables introduced in Eq. (77), w is 



the vector of classical representations of the collective 
quadratures. 



{x+,p^,x^,p+y 



(105) 



dW is a vector of Wiener increments corresponding to 
the measurement noise, 



dW ^dW {1,1,1, ly 



(106) 



and dw is a vector of independent Wiener increments 
describing the noise associated with the coupling of the 
oscillators to their mechanical environments, 



dw = {dWi^ , dWp_ ,dWi_, dWp^ f 



(107) 



Now Eq. ( 103[ ) is, in general, an eif^/it-dimensional 
Ornstein-Uhlenbeck process, with the steady-state vari- 
ances given by solutions of the Lyapunov equation. 



SH + HS^ = TT^, 



(108) 



where H denotes the steady-state covar iance matrix in 
the ordered basis Y. Solving Eq. (108) provides all of 



the conditional and unconditional variances of interest. 
Note, however, that the matrices S and T themselves 
depend on the steady-state conditional variances, which 
must first be obtained by solving the algebraic Riccati 
equation giving the steady-state of Eq. (81). 



2. The Symmetric Case 



In the fully symmetric case, the e quati ons for the mea- 



sured subsystem (X_|_,P_) in Eq. (1031 decouple from 



those for the perturbed subsystem (X_, P+). Since = 
now, the feedback force applied to both mechanical os- 
cillators, from Eq. (101), should be 



F{t) 



71 



Xj_ sina;„,i. 



(109) 



The equations for the measured subsystem alone can 



be written in the same form as Eq. (103), with the 
appropriate matrices specified in App. F3 This now 
describes a /owr-dimensional Ornstein-Uhlenbeck pro- 
cess, with steady-state solutions for the covariances again 



(103) given by the solution of Eq. (108). The reduced system 



size in this case facilitates a perturbative solution. Ex- 
panding in the (assumed) small parameter (1 -I- a)~^, 
where a is the feedback gain in Eq. (109), we find 



vt 



Vx^ + AtjC 



V- 



X. 



V^_ ^ Vp_ + 477C 



^4 



1 



(110a) 
(110b) 



to first-order in (1 + a)^^ and in the regime $1^7. That 
is, in the limit of a large feedback gain a, the uncondi- 
tional variances reduce to the conditional variances. The 
conditional two-mode squeezing and entanglement cal- 
culated in Sec. |V] can therefore be converted to uncon- 
ditional two-mode squeezing and entanglement via feed- 
back. In the opposite limit of no feedback (a — 0), we 



have y. 



fb 



yPo 



nth 



1/2, as expected. 



Note that in the (opposite) case il = 0, the result 
of Eq. (110a) is exact, consistent with the known result 
for single-mode back-action-evading measurement [6J . In 
this case we also have V^_ = Vp_ , as the observable P_ 
is neither measured nor perturbed by the measurement. 



3. The Asymmetric Case 

Again, in the asymmetric case all conditional and un- 
conditional covariances may be obtained by solving the 
system of Eq. (108). However, the size of the system is 



such that we cannot obtain useful analytical solutions, 
and we resort to numerical solutions. Fig. [7] shows un- 
conditional variances under feedback as a function of the 
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7: (Color online). The steady-state unconditional vari- 



ances under feedback, V-z^ and Vp , as a function of the feed- 



back gain a (introduced in Eq. ( 101 1), for an optomechanical 
coupling asymmetry of Gd/G = 0.05, nth = 25, G — 500, and 
in the regime n 3> 7. The steady-state conditional variances 
V-^ and Vp are shown as the dashed horizontal lines. These 
conditional variances are unequal due to the coupling asym- 
metry. The steady-state unconditional variances approach the 
steady-state conditional variances asymptotically as the feed- 
back gain is increased, though not as rapidly as in the case of 
symmetric coupling. 

feedback gain, in the uncompensated asymmetric case 
and for an optomechanical coupling asymmetry of 5%. 
It is seen that these steady-state unconditional variances 
{V^ and Vp^ ) approach the steady-state conditional 

variances {V-^ and Vp ) asymptotically as the feedback 
gain is increased. However, this approach is slower than 
in the fully symmetric case. In the compensated asym- 
metric case, the unconditional variance V^ is the same 

as Vp' , and both tend to the (common) steady-state con- 
ditional variances under strong feedback. 



D. Summary 

In this section the unconditional variances of the col- 
lective mechanical quadratures were calculated. The fact 
that there is two-mode squeezing in the conditional vari- 
ances does not imply that there is two-mode squeezing 
in the unconditional variances. However, feedback can 



be employed (see Eq. (101)) to reduce the unconditional 



variances to their corresponding conditional variances 
(see Fig. [7]). 



VII. CONCLUSIONS 

A back-action-evading measurement of a collective 
quadrature of two mechanical oscillators may be per- 
formed using a cavity detector. This is achieved, for 
a system operated in the resolved-sideband regime, by 
two-tone driving of the cavity with a detuning of plus 
and minus the average of the two mechanical oscillator 



frequencies. 

For the purpose of force sensing, one can surpass the 
full quantum limit on resonant force sensing provided 
the damping rate of the auxiliary oscillator is much lower 
than that of the driven oscillator, and the coupling asym- 
metry is matched accordingly. This would be experi- 
mentally challenging, so perhaps more useful is the fact 
that one can significantly surpass the standard quantum 
limit (and also surpass the full quantum limit) on de- 
tuned force sensing without such constraints on system 
asymmetries. 

Further, in the adiabatic limit and in the regime 
where the collective mechanical oscillator frequency is 
much greater than the average mechanical damping rate, 
one can conditionally prepare an entangled two-mode 
squeezed state of the mechanical oscillators via measure- 
ment of the cavity output field. The presence of this en- 
tanglement may be verified in a straightforward manner 
from the measurement record. Further, simple feedback 
based on the measurement record may be used to convert 
the conditional two-mode squeezing and entanglement to 
unconditional entanglement. 

Given the rapid progress in the measurement of mi- 
crowave fields [SI [Sni [57], including demonstrations of 
feedback control [5H| HHI, the proposed scheme appears 
to be a promising route towards the generation and veri- 
fication of macroscopic, all-mechanical entanglement and 
force sensing beyond conventional quantum limits. 
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Appendix A: Derivation of Hamiltonian 

We start from the Hamiltonian of Eq. ([5]), with the 
signal force given by Eq. ([m. Moving into an interaction 
picture with respect to Hq — curnia^o, + wb) + uicc'c, we 
can rewrite this Hamiltonian as 

n ^ n (a'fa - Pb\ + ga (ae~^"'"* + a^e*"''"*) c^ c 



+-ffdr 



(Al) 
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The environments of the three oscillators are assumed to 
be ensembles of non-interacting oscillators. The usual 
Born, Markov and rotating-wave approximations are 
made on the system-environment interactions, and we 
also ignore environment-induced level shifts, as is typical 
for quantum optical master equations [4]. Therefore, the 
Heisenberg equations corresponding to Eq. ( Al I , neglect- 
ing noise terms, are 



h = inh-igte^^-^^c^c- —b, 



-igb [be 



6te«"".* I c 



-if_, e"'""* - i£:_e'"'"* 



(A2a) 
(A2b) 

(A2c) 



Adopting the ansatz c{t) = co(i) -I- c+(t)e~"^'"* + 
c_(t)e"^™*, assuming w„j >• k (resolved-sideband 
regime) , explicitly separating out the Fourier components 
of the field at the driven sidebands, and then equating 
frequency components [BO], we obtain the system: 



-in 



a - iga i 



4c+ + cLcq 



la. 



b = iflb- igh [clc+ + c^Coj - yS, 
Co — —iga (ac_ + a)c+) — igt (6c_ + Pc+ 

-2C0, 
c+ = -it+ + iLOmC+ - -c+ ~ igaaco 

-igbbcQ, 
c- = -i£- - iuJmC- - -c_ - igaa^co 

-igbPcQ. 



(A3a) 
(A3b) 

(A3c) 

(A3d) 
(A3e) 



Solving Eqs. (A3d) and (A3e) for the steady-state at 



the driven sidebands, assuming the optomechanical cou- 
plings are relatively small, yields: 



(fi+> 



-184 



-iujm + '*/2 
-i£^ 



c+, 



c_ 



(A4a) 
(A4b) 



iujm + k/2 
Then we can write down the effective Hamiltonian. cor- 



responding to Eqs. (A3a)-(A3c), as 
n ^ n (a'' a - b^'b) + ffdiss 



+9a {c^c'a' + cLca' + c*^_ca + c^c'a^ 

+gb (c+c'b'^ + cLcb^ +c\cb + C-C^b\ . (A5) 

Now we assume that the cavity sideband amplitudes 
will have the same steady-state amplitudes but different 



phases, c+ = ce~*'^ and c_ = ce^^ . Here c is assumed to 
be real, without loss of generality because it only changes 
the cavity quadrature to which the collective mechani- 
cal quadrature is coupled, and is readily compensated 
by adjusting a local oscillator phase. The Hamiltonian 
becomes: 



n = n 



a^a - Pb\ 



"^QaCXa^^Xc + 2gi,cXi,^i,X, 



+^diss, (A6) 

where we have introduced quadratures of mechanical and 
electromagnetic modes as Xq^^ = {qe^^+q^e~^'^)/V2 and 
Pq_^ = —i{qe^^ — q^e~^^)/V2, with the lack of a phase 
subscript implying that the phase has been set to zero. 
The Hamiltonian of Eq. (|A6| may be rewritten as 



n = n(x+X^+P+p\ +G(cosipX+-smijP+]Xc 
-Gd (cos Vl_ - sin^P_) Xc + ^diss. (A7) 

Irrespective of the value of ip, there exists back-action 
in the "unperturbed" subspace due to asymmetric cou- 
pling {Gd 7^ 0), and so we may set V' to zero. Then the 
Hamiltonian reduces to that given in Eq. ^. 

Appendix B: Heisenberg-Langevin Equations 

The system of Heisenberg-Langevin equations describ- 
ing the dynamics of the two mechanical oscillators is 
given in Sec. In] by Eq. (12 1, repeated here for conve- 
nience, 

"^ry^M-y + FBA + N-e, (Bl) 



dt 



where V is the vector of collective mechanical quadrature 
operators of Eq. (Ill and ^ is the vector of input noise 
operators of Eq. ( 15 ). The coupling of the mechanical os- 



cillators to the cavity mode is included in Eq. (Bl ) as an 
inhomogeneity, the vector of back-action forces Fba- In 
the fully symmetric case (symmetric mechanical damping 
rates and symmetric optomechanical coupling rates) the 
system matrix M is given by Eq. (13), the back-action 
force vector -Fba is given by Eq. ( 14 ) , and the input noise 
weighting matrix N is given by Eq. (|15[) . The Heisenberg- 
Langevin equations for the cavity quadratures themselves 
are given by Eqs. (16a) and (16b). Here we give the 



Heisenberg-Langevin equations describing the system in 
the presence of both damping and coupling asymmetries. 
These equations retain the form of Eq. (Bl ) in all cases, 



and they shall be specified both in the original and ro- 
tated bases, and both without and with the compensating 
parametric driving. 



1. Original System - Original Basis 

In the presence of both coupling and damping asym- 
metries, in the basis of the original collective observables 
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of Eq. ( 11 ) and without compensation, the system matrix and the vector of input noise operators as introduced in 



M = 



-7/2 n -d-//2 

-n -7/2 -d'y/2 

-d-i/2 -7/2 VL 

-d7/2 ~VL -7/2 



(B2) 



where the average mechanical damping rate 7 and the di- 
mensionless damping asymmetry d are as introduced in 
Eqs. (pSl) and ptl), respectively. Clearly from Eq. (lB2l), 



the asymmetric damping directly couples the nominally 
"measured" (X_|.,P_) and the nominally "perturbed" 
(X_,P+) subsystems. The vector of noise input opera- 



tors is given by Eq. (151, though the noise input weighting 
matrix is 



N = 



v/^^M v/^TIdJ 







Vt^M 







(B3) 



where we have introduced the collective damping rate 
notation, 



Tt^ = V7 (VTTd ± ^r^d) /2. 



(B4) 



The vector of back-action forces associated with the cou- 
pling to the cavity is now 



Fba = (0, GdX,, 0, -GX,f 



(B5) 



with the asymmetric coupling rate Gd as introduced in 
Eq. (10b ). Th e equation for the cavity quadrature X^ is 
still Eq. (16a), while the equation for P^ is now 



Pr 



GXj^ + GdX— 



:Pr 



^Pr 



(B6) 



According to Eq. (B6), the cavity is no longer coupled to 



(i.e. measures) X+ alone, but to a linear combination of 
X+ and X_. Accordingly, the observable P_ (which is 
dynamically coupled to the measured observable X^) is 
now heated by back-action from coupling to the cavity, as 
is also clear from Ec^. ( B5 ) . Clearly, both types of asym- 



metries ruin the perfect back-action-evading property of 
the measurement scheme. The significance of these asym- 
metries is addressed throughout the paper. 



2. Original System - Rotated Basis 

In the presence of optomechanical coupling asymme- 
try it is convenient to move into a basis of rotated col- 
lective quadrature observables, as defined in Eqs. (18a) 



and (18b). The Heisenberg-Langevin equations retain 



the form of Eq. (Bl ), with V now denoting the vector of 



rotated quadrature observables as defined in Eq. (|21|) 

V 



(X+,P_,X_,P+ 



(B7) 



Eq. (22) 



S — (^+,in7 ^-,irn ^-anj ^+,in) i 



(B8) 



with the rotated collective input noise operators defined 
in terms of the original collective input noise operators 
in the obvious manner. 



The system matrix is now, c.f. Eq. (B2|, 



M = 



-7/2 h -d"f/2 -pQ. 

-n -7/2 -pfl -d'y/2 

-d-i/2 pii -7/2 n 

pVL ~d-i/2 -n -7/2 



(B9) 



with fi as introduced in Eq. (20b), while the matrix N is 
still given by Eq. ( B3 1 . The back-action force vector is 
now 



PBA = (0,0,0,-GXe 



(BIO) 



while the Heisenberg-Langevin equations for the cavity 
quadratures are now Eqs. (16a) and, c.f. Eq. (B6), 



Pc — — GrA+ — zz Pc 

^2 



V^Pc, 



(Bll) 



From Eq. (B9), it is clear that in the rotated basis both 



the damping and the coupling asymmetries directly cou- 
ple the nominally measured (X_|_ , P_ ) and nominally per- 
turbed (X_,P_|_) subsystems. However, now only the ro- 
tated observable X+ is coupled to the cavity (i.e. mea- 
sured); indeed, this motivated the definition of the ro- 
tated observables in the first place. Accordingly, and as 



seen in Eq. (BIO), the back-action of the cavity only di- 
rectly heats P4.. 



3. Compensated System - Rotated Basis 

Including compensating parametric driving as per 
Eq. (25), the system matrix in the rotated basis of 
Eq. (|2lf is 



M 



-7/2 n -di/2 

-^ -7/2 -^7/2 

-d-i/2 2pQ. -7/2 Q. 

2pVL -d'y/2 -n -7/2 



(B12) 



The effect of the compensation is that the nominally 
measured subsystem is no longer coupled to the nomi- 
nally perturbed subsystem via the coupling asymmetry, 
though it remains coupled through the damping asymme- 
try. This partial decoupling shall prove to be useful. The 
other parts of Eq. (Bl I remain unchanged: the vector of 



noise input operators is given by Eq. (22), the vector of 
back-action forces is given by Eq. (BIO), and the equa- 



tions for the cavity quadratures are given by Eqs. (16a) 
and (|B11|. 
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Appendix C: Calculation of Spectra 

1. Heisenberg-Langevin Equations in the 
Frequency Domain 

The back-action of the cavity on the mechanical oscil- 
lators may be determined by calculating noise spectra of 
the measured and perturbed observables. The system of 
Eq. (Bl I is readily solved in the frequency domain in all 



cases. Taking the Fourier transform of Eq. (Bl) leads to 



VM = -xM • N • CM - xM • ^baH, (ci) 

where we have introduced the susceptibility matrix, 

xH-(M + zwl)-i. (C2) 

For convenience, we also introduce the thermal suscepti- 
bility matrix, 



X[w] = xM ■ N. 



(C3) 



The symmetrized noise spectral densities of i nte rest, for 
an observable Z = [V]i, are defined in Eq. (26). They 



may be obtained from frequency-domain solutions of the 
Heisenberg-Langevin equations as 

(zM)tiM. 



■S'zH 



(C4) 



From Eq. ( CI ), there are clearly two contributions to the 



noise spectra of the collective mechanical quadratures. 
The first term describes the intrinsic thermal and quan- 
tum fluctuations of the mechanical oscillators, while the 
second term describes the back-action heating of the me- 
chanical oscillators due to their coupling to the cavity. 
We now calculate these contributions in turn. 



2. Thermal Noise Spectra 

The thermal noise spectra ultimately depend on the 
correlation functions of the input mechanical noise oper- 
ators; these are 

(^^in(t)^«,i„(t')> - ^{na + nb + l)S{t-t'),iC5si) 

(Z,,i„(t)4i„(t')) = ^{na-nb)6it-t'), (C5b) 

for i =/= j in the second correlation function, and where 
Z G {X,P} and i,j G {+,—}. Note that na{nb) is the 
thermal occupation of the environment to which mechan- 
ical oscillator a{b) is coupled. Therefore, the thermal con- 
tribution to the noise spectral density of the measured 
observable X+ is 



S 



■ 2 

Yl (xJm[HXl(m+2)H 

.771—1 

+Xl(m-H2)HXlmH)J , (C6) 



+ 2(«a 



n-fcj 



with Xjfe['^] denoting elements of the thermal suscepti- 
bility matrix introduced in Eq. ( C3 ) . This quantity is 
evaluated explicitly, both at the resonant peaks (Q in 
the original case, n in the compensated case) and at a 
large detuning A, in Sec. IIIB 



3. Back-Action Noise Spectra 

The contribution to the noise spectral density of the 
measured observable due to back-action is itself depen- 
dent on the dynamics of the cavity. Accordingly, we first 
take the Fourier transform of the coupled cavity quadra- 
ture in Eq. (16a), 



where we have introduced the cavity susceptibility, 

■-^[uj] = -iuj + K,/2. (C8) 



Xc 



The input noise correlation function for this cavity 
quadrature is 

(l,,in(t)^c,in(i')> = I (2nc + 1) 6{t - t'), (C9) 

where fie is the thermal occupation of the cavity bath. 
Frequently, in experiments, fic will be very close to zero. 



From Eqs. (BIO I, (CI ), (C7), and (C9), the back-action 



contribution to the spectrum of the measured observable, 
assuming that k^ fl, is 



SfU^j\Xi4M'C{2n, + l), 



(CIO) 



where X14M is an element of the susceptibility matrix 
introduced in Eq. (C2), and C is the rotated coopera- 
tivity parameter introduced in Eq. (l36|). Eq. (CIO) is 



evaluated explicitly at resonance and far-detuned from 
resonance in Sec. IIIIBI 

However, for the purpose of assessing the system's use- 
fulness for force sensing, we need to know the bandwidth 
over which quantum limits can be surpassed. This re- 
quires knowledge of the full frequency-dependent contri- 
bution to the noise spectrum. In general this is com- 
plicated, though reasonably simple forms may be ob- 
tained in the perturbative regime (i.e. to second-order 
in p and d), and in the extreme asymmetric case (d = 1 
and p = —1 without compensation, or just d = 1 in the 
case with compensation). 



4. Output Noise Spectrum 

The total output noise spectrum of the measured col- 
lective mechanical quadrature consists not only of the 
quantum and thermal fluctuations of the mechanical os- 
cillators (calculated in App. C 2 1 and the back-action 



heating due to the coupling to the cavity (calculated in 



23 



App. C3), but also includes the noise added by the detec- 
tor (i.e. the cavity). The total noise spectrum is readily 
calculated using the input-output formalism of quantum 
optics [3], and this will facilitate comparisons with con- 
ventional quantum limits on measurement [SI. 

To calculate this noise spectrum, we need to know how 
the collective mechanical quadrature couples to the cav- 
ity field. This is given, in the frequency domain, by the 
Fourier transform of Eq. (Bll I, 



facilitates a comparison with conventional quantum lim- 



p,M = -xcHGl+H + xcHT^AanH, (cii) 

where the cavity susceptibility Xci'^] is defined in 
Eq. (ICSl). From Eqs. ((C7| and (|C11|), we have for the 



cavity mode annihilation operator, 



(C12) 



Applying the usual boundary condition for a single-sided 
optical cavity. 



we find the cavity output field to be 

i{uj — LOc) + k/2 , 



(CIS) 



Cout W 



-Couti^ 



-iXc(w - 



i{io — uJc) — k/2 
uJc)\/k/2GX+[uj - Wc], (C14) 



with the frequencies now specified in the (non-rotating) 
laboratory frame. Here Cout ['-^] describes the output field 
due to the coherent cavity driving fields, the second term 
describes fluctuations of the input field filtered by the 
cavity, and the third term carries the signal due to the 
mechanical oscillation. 

It is assumed that the output cavity field is subject 
to homodyne detection in the usual manner ^. The 
measured homodyne current is then 

m = blo{t)c,nt{t)+b^o{t)cLit), (C15) 

where the local oscillator amplitude is 6lo(^) — iBe^^'^"*', 
with B assumed real without loss of generality. The 
Fourier transform of each component in Eq. (C15) is 



evaluated as a convolution integral, allowing us to calcu- 
late the spectrum of the homodyne measurement current 



Si[uj], as quoted in Eqs. (30 1 and (31) [5] 



Appendix D: Force Sensing Transfer Functions 

The ability to perform a back-action-evading measure- 
ment of an oscillating mechanical observable suggests the 
possibility of force sensing beyond conventional quantum 
limits. Given that the noise spectral density of the mea- 
sured observable has been calculated (see Sec. Ill), the 



next task is to calculate the transfer function relating the 
signal force on one mechanical oscillator to the measured 
collective mechanical quadrature. This allows one to de- 
termine the noise added by the force sensing scheme, and 



its (see Sec. IV I. These transfer functions are calculated 
here. 

Suppose that the mechanical oscillator a is subject to 
the signal force /(t), as per Eq. (47). In a frame ro- 
tating at the average mechanical frequency tUm (i.e. the 
same frame used to write the Hamiltonian in Eq. (191) 
the signal force is described by 



HF = .ma + r{t)a^ 



where we have 



m = f{t)e 



-iojat iflt 



(Dl) 



(D2) 



The Hamiltonian driving of Eq. (Dl) adds a vector of 



driving forces F{t) to the Heisenberg-Langevin equations 
of Eq. ( 12 1, given by 



F{t) 



(cos 6* -sin 61) Im /(i) 

(cos 6* + sin 61) Re f{t) 

(cos e + sin e) Im f(t) 

(cos6l-sin6l) Re f{t) 



(D3) 



where 6 is the angle describing the coupling asymmetry 
introduced in Eq. (17). In order to evaluate the force 



sensing transfer function, we can neglect the noise terms 
and the coupling to the cavity in Eq. (12), leaving the 
system of Heisenberg equations. 



dt 



V = M-V + F{t), 



(D4) 



where V is the vector of rotated quadrature observ ables 
inEq. (21 ), and the system matrix M is given by Eq. (B9) 



or (B12) for the uncompensated or compensated cases, 
respectively. The system of Eq. (D4) with the drive of 
Eq. ( |D3| is readily solved in the frequency domain, lead- 
ing to the required transfer functions. These results are 
now given, for the cases of a signal force resonant with the 
mechanical oscillator and for a signal force far- detuned 
from the mechanical resonance frequency. 



1. Detection of mechanically resonant force 

First we consider a mechanically resonant signal force. 
For the uncompensated system, assume the signal force 
f{t) is contained in a narrow bandwidth about the me- 
chanical resonance frequency uja- This implies that in 
the rotating frame f[uj] is peaked at zLil. For the com- 
pensated system we assume f{t) is contained in a nar- 
row bandwidth about uj'^, leading to f[uj] being peaked 
at ztCl. The transfer function between the signal force 
and the measured observable takes the form specified by 
Eqs. (IsTl) and (Isll). 



The transfer functions are modified in the presence of 
asymmetries by the function gr{p, d). For the uncompen- 
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sated system, this modification is given by 



9rip,d) = 



2(l-d2+p2) 

/arctanp\ 
xcosec I 



-; 



(D5) 



This factor is just 1 for vanishing asymmetries {p, d = 0), 
implying that X^ responds resonantly to the force as ex- 
pected. In the absence of coupling asymmetry, we find 
g^[p = 0, d) = 1/(1 + d). That is, the gain in Eq. (Im]) is 



given by the damping rate of the driven oscillator. Ac- 
cordingly, for d — > — 1 (the damping rate of the driven 
oscillator goes to zero), the gain diverges. For the com- 
pensated system, the transfer function modification is 
given by 



9r{p,d) 



1 



I- d^ + d^p^ 



(1 — d) cos 



/arctanp 



. /arctanp\ 
-dpsm(^^j 



(D6) 



2. Detection of mechanically non-resonant force 

The second case of interest is that in which the signal 
force is far-detuned (by an amount A 3> 7) from the 
mechanical resonance frequency. The transfer function is 
now specified by Eqs. (51 1 and (65). The modifications 



to the transfer functions due to asymmetries are given by 
the function gn(p,d). In the uncompensated case this is 



1 + p - V^l + p2 /arctanp' 
9n{p,d) = ^- cosec ( ^ ) , (D7) 



while in the compensated case we find 



9n{p-,d) 



cos 



arctanp\ 
2 ) 



(D8) 



In both cases the modulation is independent of the damp- 
ing asymmetry, and so approaches 1 for vanishing cou- 
pling asymmetry. 



in the form of Eq. ( 79 1 , but now with the matrices 



M = 



-7/2 n 
-o -7/2 



(Ela) 
(Elb) 



where T is the measurement rate introduced in Eq. ( 75 1 



2. Conditional Variances 

Also, in t he fully symmetric case the system matrix M 
of Eq. (B9) simplifies to that given in Eq. (13), and the 



system of conditional covariance equations of Eq. (81) 
becomes, writing out the equations for the indepenc 
elements explicitly. 



ent 



Vx, 


- 2r!S+_ - ^Vx^ - 


Ar^VVi^ 


+ 7"tot 


(E2a) 


Vp_ 


= -2r2E+_ - 7V>_ 
+7ntot, 


-Ar^V^l 


- 


(E2b) 


Vx^ 


= 2r!S_+ - 7Fx_ - 


- ^V^^lx 


+ 7"to 


,(E2c) 


yp. 


= -2r!E_+ - 7V^P^ 

+r, 


- AriTY.\ 


f +7"t 


Dt 

(E2d) 


S+- 


= -nvx+ + ^Vp_ - 


-7S+- 




(E2e) 


S-+ 


= -0T6f_ + VlVp^ - 
-47yrSxxS++, 


-7S-+ 




(E2f) 


S++ 


-4j/rE++^x+, 


-7S++ 




(E2g) 


i;__ 


= -n-i^xx + ^^pp 

-47yrE+_Sjfx, 


-7S__ 




(E2h) 


Sxx 


-4r/rVx^Sxx, 


l^xx + 


ind 


(E2i) 


Yjpp 


= -nY.++ - sis__ 


- jT,pp + "ffid 






-4r;rE+_I]++, 






(E2j) 



Appendix E: Best Estimates of Quadratures and 
Conditional Variances - The Symmetric Case 

1. Quadratures 

In the fully symmetric case (equal optomechanical cou- 
pling rates and mechanical damping rates), the measured 
subsystem ( A"+ , P_ ) decouples from the perturbed sub- 
system {X^,P^) in Eqs. (99) and (Bl). For the measured 



with TT-tot, accounting for quantum and thermal fluctua- 
tions, as defined in Eq. (|89|). From Eq. (|E2d|), note that 



subsystem alone, we can still write the filter, for instance. 



the observable that is perturbed by the measurement, 
i-+, is heated at the measurement rate F, as required 
by Heisenberg's uncertainty principle. Further note that 
Eqs. (E2a|, (E2b) and (E2e), including the equations for 
the variances of most interest {Vx^ and Vp_), form a 
closed system and analytical steady-state solutions are 
readily found. 

First consider the simplest case, 17 = 0. Writing the 
measurement rate in terms of the cooperativity parame- 
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ter, r = 7C, we have 



Vx^ = 



V4?7C(nth + 1/2) + 1/4 - 1/2 
4r/C 

V^th + 1/2 1 



2V^ VC' 
.+ - nth + l/2 + C, 
Vjf_ = ^P_ = nth + 1/2. 



Vp. 



(E3a) 

(E3b) 
(E3c) 



The variance of the measured observable {X^) falls with 
a stronger measurement (higher cooperativity) , while the 
conjugate observable (P+) is heated to a greater extent. 
The uncoupled observables (X_ , P_ ) are unaffected by 



the measurement. Note that if we set the thermal occu- 
pation of each mechanical oscillator to be the same, the 



results of Eqs. (E3a) and (E3b) are consistent with that 



for a back-action-evading measurement of a quadrature 
of a single mechanical oscillator ^6j . 

However, for both generating entangled states and 
force sensing beyond quantum limits, we require that 
the collective mechanical oscillation frequency, Vi, is non- 
zero. Solving for the steady-state of Eqs. ( E2a[ ), (E2b) 
and ( E2e I in this case leads to a fourth-order polynomial 
This may be solved analytically, and 
._ and Vp_ then follow in turn. We 



equation in Vx_f_ ■ 
expressions for S 
find: 



Vx^ = 



Vp_ 



^ (-1 + -^ v^87772Cfitot + 7' - 41^' + Vt" + 4r!2y72 + m"^ + l&i^^^Cn,,\ 



4ryC 

1_ 

2n 



Vx^ + ^vCVl 



"-tot 



2ntot-T6f+-47?Cy|^+ryC(^) {n,,t ~Vx^ - '^vCV^^) , 



(E4a) 
(E4b) 
(E4c) 



with ntot given by Eq. (89). The results of Eqs. (E4a)- where the appropriate matrices are 



( E4c I also apply to the compensated asymmetric system, 
with the understanding that the observables referred to 
are then the rotated observables, and we must replace pa- 
rameters without tildes by the corresponding parameters 
with tildes. More useful results based on Eqs. (E4a)- 



V E 1 Note also that the second 
) follows as a limit of Eq, 



(E4a) 



( E4c I are given in Sec. 
expression in Eq. (E3a 
in the regime fi ^ 7. Further, substituting the full re- 
sult for Vx^ from Eq. (|E4a| into the expression for Vp_ 
E^ 



in terms of Vx^ from Eq. (E4c), and re-expressing the 
result in terms of Vx^ , we find the result of Eq. ( 88 ) . 



In the compensated case (with no damping asymme- 
try), the solutions for the steady-state covariances are 
given by Eqs. (E4a)-(E4c), with the understanding that 



they are now solutions for the rotated observables, and 
parameters without tildes must be replaced by the cor- 
responding parameters with tildes. Consequently, the 



asymptotic results of Eqs. (86 1 and (87) remain valid 



Gi ven the solutions in Eqs. (E4al-(E4c), Eqs. (E2g| 
( E2j ) form a linear system 



Ri 

Ai 



(E_|_+,E — ,Sxjf,Spp) , 

-7 -n n 

--/ ~n n 

n n --f 
-n -n -7 



(E6a) 



fill = 



-4777C 






-ifid 



Vx+ 












!]+_ 
Vx^ 




B12 = Bii-2pn 



Vx, 
Vp_ 



(E6b) 



(E6c) 



(E6d) 



A, • R, 



Bi. 



where Bn or B12 is the appropriate inhomogeneity in 
the symmetric or compensated case, respectively. Subse- 



(E5) quently, Eqs. (E2c|, (E2d| and (E2f I form a hnear system 



26 



of the form of Eq. (E5), now with 



i?2 = (Vx_,yp^,I]_ + f, 



A2 = 



-7 





2n 





-7 


-2n 


-n 


f7 


^^ 



B21 



Ari-fCJ:j^x - 7"tot 
2E__ 



^22 = B2i-2pn 






(E7a) 
(E7b) 

(E7c) 

(E7d) 



again with B21 or ,622 being the appropriate inhomogene- 
ity in the symmetric or compensated case, respectively. 
Note that in the compensated case, as compared with the 
fully symmetric case, there are simply additional inhomo- 
geneous terms, leading to excess heating in the perturbed 
subsystem. 

The solutions to these linear systems take simple forms 
in the fully symmetric case, particularly in the physically 
relevant case n^ = (both oscillators at the s ame t e mper- 
ature). The linear systems defined by Eqs. (E6a)-(E6d) 



and (E7a|-(E7d) are readily solved in turn, leading to 



_,!]__, E 



and subsequently to 
Vx_ = nth 
Vp^ = fith 



XX, 

C 
C 



jpp 



2n^ 



= c 



72 + 4f]2 ' 

72 + 2rj2 

72 + 4f]2 ' 

7 n^ 



(E8) 

(E9a) 
(E9b) 
(E9c) 



In the limit 7 > fi, Vx_ ^ nth + 1/2 and Vp^ ^ 
nth + 1/2 + C. In this limit, the measurement is of X^ 
alone, and so the conjugate observable P+ is heated by 
an amount corresponding to the cooperativity. In the 
opposite (and more interesting) limit, $7 ^ 7, we find 
Vx^ ^ Vp_ ^ nth+l/2+C/2. In this case, X_|_ is directly 

monitored, but it is dynamically coupled to P_. There- 
fore, both of these collective quadratures are effectively 
measured, and the corresponding conjugate quadratures 
are heated by an equal amount. 



Appendix F: Unconditional Variances and Feedback 
1. Unconditional Variances 



in the asymmetric case, provided we consider the regime 
where fl ^ "f. In the original (uncompensated) case we 
find, assuming that the asymmetries are small. 



(X^) 2nCV^ {l+p'){l + dp)-dy2 

(A + ) - 2^CV^^ (l+p2)(l_rf2+p2) ' 



(PI 



2 l+p^-d^/2 



(Fla) 
(Fib) 



Note that there is a discrepancy between Eqs. (Fla I and 



(Fib) only when there is both a coupling asymmetry and 



a damping asymmetry. 



2. Feedback - The General Case 

The equation describing both the best estimates of 
the collective mechanical quadratures and the classical 
representations of these observables takes the form of 
Eq. (103). In the general (asymmetric) case, the mea- 



sured and perturbed subspaces are coupled, such that 
the state vector contains all collective quadratures and 
the noise vector includes measurement noise and the 
noise due to the independent mechanical environments, 
as specified by Eq. ( 104 1. The system matrix in Eq. ( 103 ) 



may be expressed in the block matrix form 



S = 



M 



F + Qi 
F 



Qi 
M 



(F2) 



where (as usual) M is given by Eq. (n3|, (B9| or (B12|, 



F is the feedback matrix defined in block matrix form as 



(F3) 



with I„ denoting the n x n identity matrix, 0„i„ repre- 
senting an m X n zero matrix, and 



" (a7/2)l2 


O22 


O22 


O22 



Qi 



(F4) 



with Q as defined in Eq. (80). Further, the input noise 
weighting matrix in Eq. ( 103 ) may be expressed in block 
matrix form as 



(F5) 



Q2 


044 


044 


V'7^totl4 



where we have introduced 

Q2 = v/4r]f Diag [Vx+ , S+_ , Sxx, S++] ■ 

3. Feedback - The Symmetric Case 



(F6) 



The expressions for the unconditional variances, in the 
absence of feedback, still take reasonably simple forms 



In the fully symmetric case we need only con- 
sider the observables in the measured subsystem, such 
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that Eq. ( 103 1 describes a four-dimensional Ornstein- 

(F7) 



Uhlenbeck process. The state vector is then 



the input noise vector is 

dW = {dW, dW, dW^^ , dWp_ ) ^ 
The system matrix is now, explicitly, 



(F8) 







7/2 + ^7/2 + 4r;rVx+ -^ -'^riTVx_ 

f7 + 4?7rE+_ 7/2 -4?7ri]+_ 

a7/2 7/2 -n 

f} 7/2 



and the input noise matrix is now, explicitly. 



(FIO) 
Of course, both Eqs. (|F9| and (IFIO ) could be expressed 



(F9) in block matrix form, as in Eqs. (F2) and (F5| 
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